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ABSTRACT 


Several problems are treated which arise in the study 
of generalized airforces in unsteady supersonic cascades 
having subsonic axial velocity. The finite cascade of 
oscillating wedges is first solved numerically. When the 
thickness is very small, the general nonlinear solution 
agrees with results based on the linear theory of charac- 
teristics, while for the single wedge oscillating near shock 
detachment it agrees closely with Carrier's exact solution. 
There is some indication that thickness effects are reduced 
by cascading. Next, the finite cascade of oscillating flat 
plates is solved analytically, to the third power of 
oscillation frequency. The generalized surface pressure is 
shown to agree with exact numerical results for moderate 
blade index, but diverges in the far field of the cascade. 
For comparison with the finite cascade, there is presented 
finally an elementary periodic solution for the infinite 
cascade. A simple relationship between the two basic cascade 
models is developed. The investigation is presently restricted 
to flow regions upstream of Mach wave reflections in the 


blade passages. 
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SYMBOLS 


speed of sound in undisturbed flow 
dimensionless pivot, measured from leading edge 
chord of blade or aerofoil 

blade spacing (fig. 3.1) 


dimensionless lateral excursion amplitude of 
bow wave 


arbitrary functions of the arguments (x-fy) or 


(22 Bm) 

arbitrary functions of the argument (&+8n) 
reduced frequency (k=wc/U) 

chordwise stagger of adjacent blades (fig. 3.1) 


chordwise stagger of adjacent leading edge 
Machelanes (fig. 3.1) 


blade index 
left running Mach wave (s=x-8y) 


dimensionless horizontal and vertical perturbation 
velocities respectively (thickness problem) 


dimensionless horizontal and vertical perturbation 
velocities respectively (unsteady angle-of~-attack 
problem) 


the quantities u, and v, immediately downstream 
of bow wave respectively 


dimensionless normal velocities immediately 
upstream and downstream of bow wave respectively 


dimensionless tangential velocities immediately 
upstream and downstream of bow wave respectively 


downwash velocity 


dimensionless normal velocity of bow shock in 
stationary coordinates 





dimensionless cartesian coordinates (x=X/c, y=Y/c) 
time (t=TU/c) 
aerofoil surface (y=Y/c) 


aerofoil pressure coefficient on upper and lower 
surface respectively 


Mach number of inlet flow relative to blading 
component of M normal to bow wave 
velocity of inlet flow relative to blading 


th blade 


angle-of-attack of n 
cotangent of Mach angle (cganien) 
orientation of bow shock wave in steady flow 


change in orientation of bow shock wave due to 
aerofoil motion 


adiabatic index 
wedge semi-vertex angle 


Pe ocmdidubLigi running Mach waves respectively 
ESE SIS (Sasa She) 


transformation variable 
State variables (thickness problem) 
defined by equ. 2.43 


dimensionless cartesian coordinates, generally 
used for nth blade with origin at leading edge 


density of undisturbed flow 
perturbation velocity potential 


potentials due to thickness distribution and 
angle-of-attack respectively 


interblade phase angle 
modified perturbation velocity potential 
total pressure loading coefficient (ACp=C," - Cy") 


oscillation frequency 
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i INTRODUCTION 


loam General Problem Area 

Significant improvements in the performance of the next 
generation of gas turbine engines seem possible with 
compressors having lighter discs and blades, with a smaller 
number of stages and increased pressure ratio perstage. 
This calls for increased tip speeds which give rise to 
Supersonic flow relative to the blades, even when the axial 
flow is subsonic. With lighter blades operating in the low 
supersonic range close to the flutter boundary, ‘the non- 
linear effects of blade thickness may also be important. 

In this regard unsteady airforce calculations are frequently 
advanced with the warning that thickness effects have not 
been considered. 

Flutter problems have thus arisen which are quite dif- 
ferent from those associated with stalled compressor blading. 
The onset of this phenomenon is characterized by a sudden 
increase in blade stress at the fundamental mode frequency. 
When flutter is encountered, it is also observed that all 
the blades oscillate at the same frequency, usually in the 
torsional mode, with some constant interblade phase angle 
ie). 

With the above features in mind, the present work aims 
to provide methods for determining oscillating airforces in 


high-speed compressor stages for use in flutter calculations. 





Both thickness effects and certain aspects of the linearized 
Supersonic flow theory are considered. In order to study 
these problems, however, the actual flow through the com- 
pressor must be simplified. Thus, irreversible shock losses 
and viscous boundary layer interactions are neglected. How- 
ever, the resulting potential flow field remains highly 
three-dimensional due to rotational effects and downwash 
induced by the blade-tip vortex system [2]. The fully 
three-dimensional case is too difficult to handle by present 
methods and further simplification is achieved by unwrapping 
the annulus of blades from the body of the compressor (see 
Big.) I.1l). Thus, generalized airforces for flutter calcula- 
tion purposes are at present calculated by determining the 
irrotational flow of a perfect gas through cscillating 
two-dimensional cascades. 

Supersonic cascades are of the two types illustrated in 
fig. 1.1. When the axial velocity into the compressor stage 
1s supersonic, the Mach lines of the effective inlet flow 
are swept back beyond the leading edge locus. This is the 
rectilinear cascade with the supersonic leading-edge locus 
condition, in which no disturbances exist upstream of the 
blade passages. The entire flow in this case is a periodic 
extension of the flow between the first two blades. The 
solution of this problem, using a Laplace transformation, 
was first given by Lane [3]. When the axial velocity is 


Subsonic there results the subsonic leading edge locus 
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ae ae vee E The two basic rectilinear cascades 


condition, in which disturbances exist infinitely far 
upstream. This is the case of main practical interest and 
it has yet to be solved in closed form for general oscilla- 


tion frequencies. 


ey 





ez PeGellineamecascade Theory - Recent Progress 


Noteworthy of the progress which has recently been made 
with the above problem are the approaches of Platzer and 
Chalkley [4], Brix and Platzer [5], Verdon [6] and Nagashima 
and Whitehead [8]. Although somewhat different procedures 
were used in each case, namely characteristics [4,5], finite 
differences [6] and Singularity distribution methods [8], the 
agreement between the methods is generally good. The 
dependent variable in these methods is the forced steady 
state amplitude of the velocity potential. However, the 
memnods Cxhibit several disadvantages. First, convergence is 
slow and it is usually necessary to add the effect of a large 
number of aerofoils; Verdon shows regular oscillations in the 
amplitude of the unsteady lift and pitching moment out to 
blade index 30. Second, for large blade index, significant 
numerical errors are introduced in pressure loading, particu- 
larly when the blades oscillate in phase. Third, no 
mathematical proof exists to show that the limiting pressures 
are indeed representative of the infinite array. 

An alternative numerical approach by Sisto and Ni_ [7] 
attacks the infinite cascade by directly introducing the 
boundary condition which ensures periodicity of the flow with 
blade index. Instead of solving the small perturbation wave 
equation for the forced steady state amplitude of the velocity 
potential, Sisto solves the complete set of unsteady flow 


equations using the time-marching technique. However, in 


12 





the example given, 1000 time-steps were needed for 
convergence [7]. The method, apparently, may be extended 
to include thickness effects. 

The infinite cascade problem has recently been solved 
analytically to the first order in frequency by Kurosaka [9]. 
Using the streamwise coordinate x as the Laplace transforma- 
euenm Variable, a general solution in integral form for (x,y) 
is obtained, which contains arbitrary functions for the 
mittial conditions ¢(0,y) and 3¢/dx(o,y). The complexity of 
the initial value problem is reduced by representing $¢$(x,y) 
as the first two terms in a power series based on frequency. 
Reformulation as a simple integral equation follows by 
eliminating the unknown initial condition using the periodic- 
ity requirement in the cascade. The analysis is complicated 
and the equations of great length. The stability boundaries, 
for torsional oscillations at reduced frequency 0.1, are of 
little practical application in turbine flutter calculations. 
Comparison with any of the available finite cascade solutions 
would have been of more fundamental interest. A limitation 
of the theory is that when the blades oscillate in phase, the 
Surface pressures become infinite; the lifting pressure, 
however, remains finite. This could indicate that a first- 
Order power series in frequency is inadequate. Kurosaka 
is presently extending this work to include terms of higher 
order in frequency.” In a very recent paper, Verdon and 


&s 
Unpublished communication. 
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McCune present a singularity method for the infinite 
cascade [10]; further work, however, iS required to clarify 
the numerical convergence aspects of the method. 

The concern of the discussion to the present juncture 
has been with linearized supersonic flow theory. Thus, the 
above methods for predicting flutter in compressor stages 
have employed unsteady supersonic airforces calculated on 
the basis of flat-plate theory; the influence of blade 
thickness remains to be investigated. Thickness influences 
the steady flow about supersonic wings by causing a redistri- 
bution of pressure which usually changes the total lift force 
slightly and the pitching moment considerably. This well 
known experimental result is in agreement with Busemann's 
second-order potential flow theory [ll]. The unsteady 
counterpart of the Busemann theory is the second-order 
potential flow theory of Van Dyke [12]. Van Dyke gives a 
solution for the slowly oscillating aerofoil of arbitrary 
profile expressed as a third-power series based on frequency, 
and a solution for the oscillating wedge which is exact in 
frequency. Teipel considers the unsteady supersonic aerofoil 
with arbitrary profile, attacking the problem numerically 
[13,14]. This method also assumes potential flow, with no 
entropy loss through the oscillating bow shock wave. However, 
shock losses may be of practical interest in high-speed 
compressor blading near shock detachment. These effects may 
fortunately be estimated by comparison with Carrier's exact 


solution for the oscillating wedge [15]. 
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i. 3 Present Investigation 

The main aims of the present work were twofold: first 
it was desired to determine a numerical procedure for 
calculating the nonlinear effects of blade thickness in 
finite cascades, and second it was desired to investigate 
theoretically the features of the linear (flat plate) 
cascade. The nonlinear study employs a generalization of 
Teipel's method, briefly mentioned earlier in this sequel. 
However, the flow variables immediately behind the second 
oscillating bow shock wave, and indeed behind each succeeding 
wave, can no longer be determined simply by assuming steady 
upstream flow conditions. The flow is taken to be irrota- 
tional and isentropic and to satisfy the unsteady form of 
the nonlinear transonic small perturbation equation at all 
points in the field. Teipel's approach is to derive the 
unsteady boundary conditions behind the oscillating bow 
shock wave (an interesting strategem involving the Rankine 
Hugoniot equation is used) and to then develop the solution 
downstream uSing a nonlinear characteristics procedure. A 
concise account is given of Teipel's work, omitting considerable 
detail, but emphasizing fundamentals and important computational 
steps. 

Before extending the basic method to the cascade problem, 
its accuracy in the low supersonic range near shock detach- 
ment was carefully investigated by comparison with Carrier's 


exact solution for the oscillating wedge. This elegant 


aD 





solution involves evaluation of an infinite Bessel function 
series; perhaps as a consequence there are no published 
results on pressure distribution, or indeed results of any 
kind showing the effect of thickness on the main aero- 
dynamic stability derivatives in the low Supersonic range. 
Such results, given here for the oscillating wedge, are thus 
of fundamental interest, as well as providing the necessary 
base for comparison with the potential flow theory. 

Agreement is found with Van Dyke's result, that thick- 
ness generally reduces the dynamic stability of the torsion 
mode. However, the magnitude of the effect on pressure 
distribution and various unsteady force eae eves! is more 
Gramatically illustrated by the present results than by those 
of previous authors. In the case, for example, of slow 
oscillations about a forward pivot at low supersonic Mach 
number, it is shown that while the biconvex aerofoil 
exhibits a level of instability only slightly higher than 
that predicted by the classical linear theory of Garrick and 
Rubinow [16], the equivalent single wedge is drastically 
unstable. 

An important limitation in scope was imposed throughout 
the present work by restricting attention to the flow in the 
cascade upstream of any Mach wave reflections in the blade 
passages. These flow regions are indicated in fig. 1.2 as 
the preinterference zone A and first blade passage zone B 


of the cascade; they frequently extend through most of the 
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fig. 1.2 Preinterference and blade passage zones 
considered in the investigation 


‘cascade. It is nevertheless to be emphasized that under no 
circumstances has attention been given to the trailing edge 
wake, which clearly influences the cascade exit flow [9,17]. 
Results for the finite cascade of oscillating wedges 
reduce correctly for vanishing blade thickness to those of 


Platzer and Verdon, mentioned earlier in this sequel. * 


* 

SUChMeG act eMInteing results, Obtained from the general 
Poumeincar thickness model, are extensively used in this 
dissertation for comparison with linear theory developments. 


IL 





Also for this limiting case, certain simple forms of the 
general bow shock wave analysis are presented. Pressure 
Gistributions are given for two 5-percent thick wedges 
peeiilating in phase for several oscillation frequencies. 
However, until the present solution is extended to the exit 
plane of the cascade the results are of limited use in 
practical flutter calculations. Nevertheless, reduced 
influence due to thickness for the second blade is indicated. 
Section 3 presents a theoretical approach, using Laplace 
transformation, to the finite flat plate cascade. There is 
developed, to the third power in frequency, an analytical 
solution for the velocity potential together with closed- 
form expressions for unsteady pressure distribution. For 
two blades, the out-of-phase (damping) solution is shown to 
be accurate for reduced frequency 0.4 at M= 1.25 and for 
reduced frequency 0.7 at M = 1.6, well beyond the useful 
range of the first-order solution obtained by neglecting the 
cubic frequency term. On the other hand, the complete third- 
order theory is satisfactory out to blade index 9 for 
oscillations with reduced frequency 0.2 at M= 1.6. The 
theory also correctly predicts, when the blades oscillate in 
phase, a strong sinusoidal modulation of pressure with blade 
index. However, owing to an expansion of the Bessel function 
which assumes small values of the argument, the third-order 
solution diverges in the far field, failing to approach a 
eenite Limit with increasing blade index. When the higher 


order terms are removed, the first-order frequency solution 
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solution is relatively simple. In paragraph 3.4 there is 
given an alternative derivation of this important first- 
order solution for the finite cascade, based on Sauer's 
classical treatment of the oscillating flat plate in unbounded 
mow (18). 

The coefficients of upper surface pressure and total 
pressure loading, according to the first-order theory, are 


found to be 


n-l 
oe -~ c() , , 2km Y . @ PH 
n L 8 xT 
eile 2 Ale Rea 4 
Ne = 2C [e L 1j] + i-— [-——-m]j e u 
P P B B 
n 1 
th 


for the n blade in the torsion mode when angle-of-attack 


is Maximum. Here C corresponds to the first blade; , is 


(u) 
ey 
interblade phase angle; k reduced frequency; d interblade 
spacing; and m the chordwise stagger of adjacent leading edge 
Mach lines (for the sonic leading-edge locus, m= 0, and 

td is independent of blade index). To first-order frequency, 
the pressure loading for the finite cascade is therefore 
constant for all blades except the first, while for p = 0 it 
is independent of position along the chord and pivot axis 


location. The upper surface pressure, however, exhibits a 


continuing oscillation with blade index, except in the case 


1M) 





of zero interblade phase angle, when the out-of-phase 
component increases by the constant amount 2km/8> from 
blade-to-blade. 

There is derived in Section 4, using a much simpler 
approach than Kurosaka [9], an elementary periodic solution 
for the infinite cascade. When the condition of flow 
periodicity is imposed, the pressure loading coefficient for 
the infinite cascade is found to agree with AC, given above, 
while for finite interblade phase angle the Beate solution 


for the surface pressure coefficient is precisely the mean 


(u) 
P e 
nN 


value of the oscillatory series for C 


1.4 Important Note 

Concerning the simple harmonic motion assumed for aerofoil 
angle of attack, the following should be emphasized. First, 
the time dependent angle of attack of the isolated blade, 
, = exp ikt, where k is reduced frequency and t nondimensional 
time, assumes an oscillation amplitude of unity. This is 
consistent with the basic premise that all flow variables are 
linear functions of angle of attack and that the onset of 
flutter is characterized by oscillations of very small 
amplitude. Where significant nonlinearities are present 
(Section 2), they enter on account of aerofoil thickness. 


th 


Second, the angle of attack of the n blade, at time zero, 


is denoted a, = O,€xP i(n-l)u , where uw is interblade phase 
angle. Solutions for surface pressure consist of a real and 


an imaginary part, and throughout this dissertation final 
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results for the ae blade are expressed at the instant when 


= ie oOyemultielyines by exp -1i(n-1)u. Accordingly, the 
real component of surface pressu.e is in-phase with the 
angle of attack of the ao blade, while the imaginary 
component is the pressure which occurs at the instant when 
the blade presents zero angle of attack to the free stream 
with its trailing edge moving down. 

This work was made possible for the author through the 


guidance of Professor Max F. Platzer of the Department of 


Meronautics. 
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2 NOME AR Tan LCKNESS EFFECTS 


The methods which are being increaSingly applied to the 
study of flutter in high-speed compressor stages employ 
generalized air forces calculated on the basis of linearized 
Supersonic flow theory. Frequently, these methods are 
advanced with the warning that the results do not include 
nonlinear thickness effects. Accordingly, there is given 
in the present section, employing a generalization of Teipel's 
method for the single oscillating aerofoil, a numerical 
solution of the thickness problem for aerofoils oscillating 
fn cascade. The method is intended for application in the 
low supersonic range near shock detachment, where thickness 
effects may be of considerable practical interest. 

Leal Aerofoil in Unbounded Supersonic Flow 
The Teipel Approach 
The small perturbation equation describing the unsteady 


transonic flow of an irrotational fluid is given by Landahl 





ier | 
2 2 Z 2 
pe -20Ggy ge] 2g ag. Fee ee ay 
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where the nonlinear term involving the adiabatic index is 
shown in ref. 19 to be of vanishing significance for 


unsteady flows in which 


Ze 





bs (Ih ea 


where k denotes reduced frequency and M,. local Mach number. 
The nonlinearity also vanishes for flows about bodies of 


small thickness, such that 


k >> 5273 ZZ 


where 6 is the larger of the quantities thickness-to-chord 
ratio or oscillation amplitude. The present aim is to 
investigate the effect of thickness on the unsteady flow in 
low-supersonic compressor stages for blade vibration 
frequencies within the critical flutter range. Previous 
work, based on linearized supersonic flow theory, indicates 
that this range extends to k = 1.5 [20]. 

A clear understanding of Teipel's treatment of equ. 2.1 
for the case of the biconvex aerofoil oscillating in unbounded 
Supersonic flow is essential before considering the cascade 
problem. Accordingly, there is presented below a concise 
account of the Teipel method. This account, based on 
translations of the German by the present writer, strongly 
emphasizes the fundamentals of Teipel's analysis, highlights 
important computational steps, but omits considerable detail. 
Further, while Teipel employs the chord length c = 2 with 
origin at midchord and works in dimensional variables, these 


conventions are not adopted in the present study. 
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Oe ae Characteristic equations 


Y 
Head-Wove in 
Steady Flow 
Head-Wove in 
Unsteady Flow 
M 
eee , 
Pineal 


| ¥(X,T) 
x 


<a MW SL enn 


— 


C 


fig. 2.1 Biconvex aerofoil in torsion mode 


With reference to fig. 2.1, the basic problems 
involved in solving equ. 2.1 are threefold. First, a solution 
must be found which satisfies the differential equation at all 
Meants within the disturbed flow field. Second, it must 
satisfy the usual linearized form of the unsteady boundary 
condition, which requires tangential flow along the aerofoil 


surface 
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Third, the appropriate boundary conditions must be satisfied 
along the oscillating head wave which separates the disturbed 
from the undisturbed flow. The remainder of this paragraph 
outlines the Teipel approach for solving equ. 2.1 subject’ to 
the tangency condition of equ. 2.3. Paragraph 2.1.2 gives 
the head-wave analysis and paragraph 2.1.3 the computational 
steps. 


Upon introducing the nondimensional variables 


= aa ee) =e 
ST et Sf eae ie co! p CU 2.4 
egqus. 2.1 and 2.3 become 
2 2 a0 Co den Doe D526 
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mere Y = Y(x,t), as shown in fig. 2.1, and y = Y/c 
Consider torsional oscillations a = exp(iwT) about 


the pivot x = b, and define a reduced frequency 


Then there may be written for the total potential 9¢(x,y,t) 
the sum of a steady thickness term b, (%ry) and an unsteady 


angle-of-attack term 9, (x,y) 


Zo 





o(x,y,t) = d, (x,y) + 6, (xryyeT Kt TS | 


whereupon the differential equation and tangency condition 





become 
2 2 
0 0 > a> 
fe? ae M2 cyan i. | — - —" = 0 2.8a 
OX Jy 
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where 7 denotes the amplitude of the forced steady-state 
velocity potential and y(x) the aerofoil profile, nondimen- 
Sionalized through division by the chord. In deriving equ. 


Pea, PI xx is neglected in comparison with 9 in agreement 


Oxx’! 
with the assumption of very small amplitude oscillations. 
The hyperbolic differential equation (2.8a) is nonlinear, 


due to the term involving the adiabatic index . However, 


the solution procedure using nonlinear characteristics 1s 
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straightforward. The potential o, (x+y) represents the 
steady-state flow (assumed to be entirely supersonic) about 
a thin aerofoil (assumed to be symmetrical) at zero angle 
Sn attack. 

The hyperbolic differential equation (2.9a) is 
linear — since oe is known. Furthermore, since the coeffi- 
cients of the second-order derivatives are those in equation 
2.8a, the characteristic directions in both problems are 
identical. The potential >, (x,y) represents the unsteady 
flow due to angle-of-attack. 

The basic approach, using the method oF characeeric— 
tics 1s now apparent. First, equation 2.8 is solved. The 
compatibility relations for 9% are extremely simple. How- 
ever, they must be applied along characteristic lines in 
the (x,y) plane with varying slope. Second, employing 
appropriate compatibility relations for Oy: the character- 
1stic mesh derived for the steady problem (in the region 
formed between the head wave and the upper surface of the 
aerofoil) is used to develop the unsteady solution. 


The characteristic lines of equ. 2.8a have the slope 


yA 


along which the state variables change according to the 


compatibility relations 


A} 





where the constant for the left running a characteristics 


(upper sign) will be denoted by c“ and that for the renlrele\ie 
B 


running 8 characteristics (lower sign) by C", as shown in 
fig. 2.2. The state variables have the meaning 
2 2 


M - 1+ (y + 1) M Us Ceca 


» 
I! 
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u = 






ead-—Wave 








NaS 


fig. 2.2 The characteristic mesh 
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where u, = 3$,/dx and v 


a = I> /9Y are nondimensional 


O 
horizontal and vertical perturbation velocities, respectively. 
Along the profile yw is known, from the tangency condition 
imposed by equ. 2.8b. The curvature of the head wave and 

the structure of the entire simple wave flow thus follow 
meemelittie difficulty from equ. 2.10 to 2.12. Observe, for 
example, that the compatibility relations demand A and u 
constant along all a characteristics, which are therefore 
straight lines in the (x,y) plane. This important result 
follows if the 8 compatibility equation is assumed to remain 


# 
valid through the bow shock wave. For then there may be 


Pietiten, referring to fig. 2.2 


By 2 Pe 57 2 
he + i a5 
3/2 rs 
Ay ts is = A 5 
3 7 2 
Ne Ha do Mb 
mere AY Occurs in the undisturbed flow and uw. = 0. There 
follows 
ds 7 dy 
Ee = es 


* e 
Oswatitsch considers the problem, and by comparison 
With exact shock polar results shows that the use of equ. 2.11 
leads to a maximum possible error in X of 6% [21]. 
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and from equ. 2.10 the result that the a characteristics 

are straight lines. The computational procedure for dealing 
with the %6 problem is now straightforward and will be given 
in paragraph 2.1.3, as mentioned previously. Using equ. 2.12a, 
the differential equation describing the by problem may be 


weLtten 


2 
0. 6 0 
ee OO be a2kN’ ——~=0 2.13 
2 2 OX OX i ~P,4 
ox oy 
whereupon, along the characteristic lines previously defined 


in equ. 2.10, the following compatibility relations can be 


shown to hold 





where u, = 9>4/9,1 Vv, = 8>,/0y are respectively nondimensional 
horizontal and vertical perturbation velocities for the 
unsteady problem. The state variable X is Known everywhere. 
In equ. 2.14 the second-order term in frequency somewhat 
complicates a numerical solution Since, as is clear, the 
potential ot itself must be determined throughout the field 
(by integrating the velocities). 

Teipel's rather involved treatment of dA/dx can be 
feopiified using a geometrical interpretation. At the 


bem O(X,y) in fig. 2.3 
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a — Characts. 





mid. 2.3 Geometrical interpretation of 3i/dx 


Xx -xX Cc a 
Cc 


However, from a result of paragraph 2.1.1, the state 
variable i is constant along an a characteristic. 


meecordingly, writing ha = X. in the above equation and using 


Cc 


an obvious approximation for the denominator 


ioe. 
y =F yee mit xg+x, 
which is simply 
ar. 1 [aa 
22 fa) ba 


Sl 





The compatibility equations for the unsteady problem thus 








become 
av u 2 2.2 
1 1 OA 2kM kM a 
du, ia + 77 a Gx te T 1 ax ——— >, dx = 10) 
2.16a 
av u 2 ae 
L 2kM kM = 


along a and 8 characteristics, respectively. Also, since 


do, = 36,/axdx + 36,/ayay 
b e 
ela + ae 2.17 
1 1 Ll Ll : 
a 


Further details of important computational steps are given 
maeparagraph 2.1.3. 
2.1.2 Bow Shock Wave Analysis 

There are two important boundaries associated with 
the , (X,Y) problem. One is the aerofoil surface along 
which the flow must be tangent and along which V1 (%,0) is 
therefore given by equ. 2.9b. The other is the oscillating 
shock wave which emanates from the leading edge of the aero- 
foil. Fluid particles in the undisturbed flow, passing 
through this oblique moving shock, will have their properties 
suddenly changed. Accounting for entropy losses, this 
change could be calculated using the shock polar [22]. How- 


ever, the flow is assumed to be isentropic and no account is 


a 





taken of irreversible shock losses. Accordingly, the flow 
variables immediately behind the head wave are derived using 
the isentropic Rankine Hugoniot relations. The more exact 
shock polar is not employed. 

Teipel's analysis of the unsteady head wave is in 
fact surprisingly straightforward, when the equation for 
the change in gas speed through a moving normal shock is 
given. Such a weak shock is shown in fig. 2.4, moving from 
left-to-right with the speed W relative to the undisturbed 
inlet flow. Denoting the velocity of the latter Wore the 
shock thus has the velocity (W + w,,) in stationary coordinates. 
The one-dimensional continuity and momentum equations then 
lead to the expression given by (2.18) for Wey the gas speed 
following passage through the shock. Here it is to be assumed 
that the velocities Wye Ww and W have been normalized through 


Givision by U. An outline of Teipel's elegant use of this 


result follows: 





W+ Wp 
Undisturbed 
flow into : 
“~~ 
shock Wy, Wo 
—_—S 





a 2, 2 a 
Wn= Wn+ 2W peau ee 


i1 Gwen Flow through normal shock 
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The locus of the oscillating head wave, as it might 
appear with the aerofoil at small positive incidence in a 
low frequency Bore lenmeal modeyeis Shown in fig. 2.5. Fora 
wedge, the orientation Nes of the head wave in steady flow 
1s constant. However, for the biconvex section shown, 
Me eee Consider P(x,y) on the oscillating head wave 
in the immediate neighbourhood of the leading edge. Then, 
notwithstanding the variation of ve with y, the approximate 


horizontal coordinate of P is 


Ie 







x = COr + 2 2 ° 19 
np” Yp Thy veche'e, 
y Head wave in steady 
flow NB all velocities 
V, normalized thro- 
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/ x 
el Iw Normal shock speed 
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Yo 
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Os Aras Geometry of oscillating bow wave 
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where the nondimensional, complex, shock excursion amplitude 
g(y) has been multiplied by exp(ikt) in accordance with the 
basic assumption involving harmonic time dependence. Omitting 


this time dependence, the slope at P(x,y) is 
—— = COE ys + eo cot ys) 2.20 
p O dy O ° 


where fig. 2.5 shows the small angle y'. Employing an 
elementary trigonometrical expansion and taking tan y' = y', 
the change in the orientation of the head wave relative to 
its position in steady flow may thus be approximated by 

y' = - sin y, se Del 
The nondimensionless velocity assumed for P(x,y) 


in stationary coordinates is 


ax 


aw 
at ikg(y) 222 


meomeeoqu. 2.19, Omitting exp(ikt). Multiplying by sin Yeu 
the normal velocity of the head wave in stationary coordinates 


is 


v' = ik sin ate g(y) ee 


Sie. 





which is to be associated directly with the shock speed 
(W + Ww) irises Etatnererore follows from equ. 2.18 
that the gas speed W, immediately behind and normal to the 


head wave iS given by 


A _ _ al ya 
Wh eo acl il , 5 2224 
(yi = Ww.) 
n 
But 
ws sinly, + y') = sin Youn y'cos Ye 2122 
whereupon 
ioe ' =! = 
-. _ 2(v'-w) A 332 ne cosy. -v') i) 
Ban yt+1 Cee Sin y 
(y+1)U sin ees O 


ZG 


Employing the Binomial theorem, and after some simplifica- 
tion, the main result of the Teipel head wave analysis may 


be written 





hs 1 2/m_* 2/m_? 
— | 
YT Trisiny, lhe = E SeReEE 5 dl oll Y 
2 IL 
+ y+L A 5 V ZL 
n 
Ui sin y 
i= 2.28 
n a 
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where the second-order terms y' and v' are given above in 
terms of the head-wave displacement g(y). The tangential 
velocity of the fluid is unchanged on passing through the 


head wave and consequently 


It 1s now required to find the dimensionless pertur- 
bation velocities 39),/ax and 8¢,/dy immediately behind the 
head wave, these being the desired boundary value quantities 
for the solution of the oT problem. Formally introducing 
the overhat notation to denote quantities immediately behind 


the bow wave 


A 


inyt+ 
Uy) (w,Siny W 


cosy) = (w Siny tw, cosy) oo Ua 


(-w,,cosy+w 


<g> 
tt 


,Siny) = (-w, cosy tw, siny,) Zs UD 
where the two groups of terms on the right of equ. 2.30 
represent fluid velocities behind the head wave (a) in 
unsteady flow and (b) in steady flow. Of course, y = Ven. - 
The remaining work is routine. Neglecting squares and 


products of y' and v', there is obtained finally 


dg ; 
cs, + 1Imog Za 


SV) 





ae a : 
Vy ny at ing Dee ANS, 
m, = — Samy 2 sine 
Al vaedl Yo Yo 
_ 2k 2 Py 
mo = y+L (1+1/M_ )siTl Vs 
232 
_ =2 2 nee 
n, = avid (cos 2Y. + 1/M, es ny Le 
_ =2k 2 ee 
no = ail (1+1/M_ ne Ot Yo Sin” Y, 


The more formidable looking equations of ref. 14 result 
when the orientation of the head wave is expressed in the 


form 


an important observation, which follows from the slope of 
the left-running characteristics in equ. 2.10, and is equi- 
valent to averaging the Mach angles before and after 
leading-edge compression. 

Equations 2.3la and 2.3lb are simultaneous ordinary 
differential equations for the unsteady head-shock profile. 
They are expressed in terms of the unknown dimensionless 
perturbation velocities Uy and Vy immediately behind the 


shock. The coefficients m and n are functions of y in general 


Shs 





because of the curvature of the steady shock boundary. 
They are known from the solution of the 6 problem. 

fal. o COMputational Procedure 

Certain mathematical aspects of the unsteady transonic 

flow about a thin aerofoil have been discussed. It is now 
desired to translate these considerations into a procedure 
for solving the problem numerically. With regard to the 
steady state or > problem, the following sequence of steps 


may be used: 


(1) Guess the point Es 





mud, 2,6 Construction procedure for characteristic lines 
and steady-state bow wave 
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(ii) There may be written, corresponding to 


mmerad and 6 Characteristics which meet at P. 


a) Wie at Sa) 
AT Wy ra Wee es 
r 3/2 Sn) ee, 3/2 + UW = oF 


from equ. 2.11, and where, as noted previously, the overhat 
denotes a quantity immediately behind the head wave. But 


from the steady tangency condition of equ. 2.12b 


eS 20 Oy 
vo 7-5 (yar) eM Te 
a 
Further, HW, = Ho and hy -- es Since the two states are on 


the same a characteristic. Accordingly 


(Mn emcomocrmiict thc a characteristic through 


ES with slope from equ. 2.10 





ole 
lI 
‘ae 
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(iv) Construct the initial segment of the 
head wave between the leading edge point aa and Py with 


slope 


which 1s equivalent to averaging the Mach lines before and 
after leading edge compression. This is not the segment of 
an a characteristic. Hence, Pi on the head wave 1s 
determined. 

(v) The state variable HW, Must satisfy the 
tangency condition and is a known function of X57 as in 
step (11). Along the 8 characteristic through P. 


+ Uy = cP 


Yy 


(Sines Ue ne 
ee 

2 in af * == uf 2 
and hence Xo and ho follow by a Simple iteration procedure. 


(vi) Guess A (thereby determining P,) and 


b 


repeat steps (11)-to-(v). 


Determination of the general mesh point Ey is obvious. 
The above sequence of steps is repeated until the entire 


characteristic mesh is formed and the 6 problem solved. 
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Some experience is useful in determining a suitable initial 
point Po and the x-wise increment A. in the case of a 
Simple wedge or flat plate aerofoil, the head wave is 
straight and the state variables constant throughout the 
flow field; consequently the complexity of the 6 problem 
1s greatly reduced. 

The same characteristic mesh applies to the unsteady 
>) problem. Here, it will be recalled, there is available 
the following information: First, the unsteady tangency 
memadition of equ. 2.9b for V1 (X,0) 3 second, the compatibility 
relations of equs. 2.l6a and 2.16b, which identify small 
p@eferential changes in the state variables u, (x,y) and 
Vv, (x,y) along the a and g characteristics, respectively; 
and third, the first-order ordinary differential equations 
(2.3la and 2.31b) for g(y), the amplitude of the oscillating 
head wave, measured from steady state. For a typical problem 
the shaded area of fig. 2.7 is of infinitesimal dimension, as 
should be clear from the information Miele LOOLnOte ~*~. 
Accordingly, the leading edge values of a) and v. are assumed 


il 


to apply throughout this small region. Let wee) c or Vv) 
for convenience, denote the dimensionless vertical perturba- 
tion velocity immediately behind the head shock at the 


leading edge. Then from equ. 2.9b 


* 
imobeatevalues fOr M= 1.2, 1 = 0.025 and M= 1.15, 


T = 0.0125 were x. = 0.0017, A = 0.0055 and as 0.0088, 
A = 0.00425, respectively; with fourteen mesh points 
Po? Por Pe eee Pat 
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fee. 2s 7 


where b is the 


Since g(o) = 0. 


(uj, tu) /2 


(vj tv.) /2 


while from equ. 
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Initial step for the unsteady problem 


ds 1/18) oe ahh jy) 


dimensionless pivot. From equ. 2.31 


Also from equ. 2.31 
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where 7 would be taken as zero for the single aerofoil in 
unbounded flow. The three algebraic equations (**) deter- 
mine Uys Vie and Sy at the mesh point Pye By similarly 
applying these equations to the points Po, Pes Pio (assume 
uy (Py) and Tees) apply at the intersection of the a 
characteristic through one with the g@ characteristic through 
Pi), the unsteady boundary conditions along the entire head 
wave may be determined. Subsequently, the compatibility 
relations (equs. 2.16a and 2.16b) give the entire unsteady 


flow field. 


2.2 Discussion of Results for Flat Plate, Wedge and 
Biconvex Aerofoils 


The linearized pressure coefficient on the upper surface 


of the aerofoil is 


Uu © 
Pp ~p 


2 


a = E i 4 = -2[u.te™**(u,+ikg,) ] Dee 
1/2pu 


where ue and u, are dimensionless horizontal perturbation 
velocities due to aerofoil thickness and angle-of-attack, 
respectively. The thickness or steady state component of 
pressure is entirely real and of no direct interest. The 
unsteady component, which varies in accordance with the 


assumed harmonic time dependence, is complex and will be 


written, omitting the term exp(ikt) 


* 
Again it is to be noted that u v, and Jy should in 


el. 
fact be denoted u, (Py), v1 (Py) and g(P,). 
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ex! = -2(u, + ik) 2.35 


The present aim is to study the effect of thickness on 
oe Therefore, it 1S appropriate to first discuss the 
meaning of the real and imaginary parts. The real part, 
pares”), is the pressure which occurs as the oscillation 
starts. It is the component of unsteady aerodynamic pressure 
which is in-phase with the angle-of-attack (aerodynamic 
stiffness). The imaginary part, reese i is the pressure 
at the instant the aerofoil presents zero angle-of-attack 
to the free stream with the trailing edge moving down. It 
is due to aerofoil rotation about the pivot and therefore 
out-of-phase with the angle-of-attack (aerodynamic damping). 

A mathematical model describing Teipel's method has been 
formulated and results compared (a) with Teipel calculations 
for a family of biconvex aerofoils in the low supersonic 
range [13], and (b) with Carrier's exact solution for the 
oscillating wedge [15]. Some disagreement was found with 
Teipel, which could not be resolved. As shown in fig. 2.8 
fOr a 2% thick biconvex aerofoil oscillating in the torsion 
mode about midchord this disagreement is in ie (kena le the 
component of pressure which is out-of-phase with the angle- 
Oof-attack. Similar results were obtained for most of the 
torsion and plunge modes considered in ref. 13. However, in 
the case of vanishing aerofoil thickness, the present model 


agrees with Teipel, and indeed with the classical linear 


45 





3-0 i \ —-— REFERENCE 11. 
\ —— PRESENT TEIPEL MODEL 
6 \ 


(u) \ out of se 
rn oo) ie 





mug, 2.8 Pressure distribution for 2% thick biconvex 
aerofoil pivoted at mid-chord 
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theory of Garrick and Rubinow [16]. The structure of 
Teipel's nonlinear characteristics mesh is based on 13-points 
across the chord. A similar number was found to be satis- 
factory for all of the unsteady aerofoil motions considered 
in the present investigation. 

As already mentioned Teipel assumes isentropic flow with 
no entropy loss through the head shock. In view of the 
envisaged application of the basic method to compressor 
blade rows in the low-Supersonic range, where irreversible 
shock losses could be important, this assumption requires 
investigation. Carrier's solution for the wedge, oscillating 
about its apex in unbounded supersonic flow, is exact in 
thickness and therefore accounts correctly for entropy losses 
in the fluid provided the amplitude of oscillation is small 
(the solution is linear in angle-of-attack). It provides 
an obvious basis for comparison with the potential-flow 
theory. 

The Carrier analysis is difficult and the final solution 
1s expressed as an infinite series of Bessel functions ae 
the first kind with certain coefficients which are to be 
determined successively from three simultaneous algebraic 
recurrence relations. Consequently, the computation involved 
1s considerable, which may account for the fact that no 
pressure distributions based on the theory are available. 
Reference should be made, however, to exact calculations by 
Van Dyke of the stability derivative my fon a 56 thick 
oscillating wedge. These are published a met. 623, which 


extends Carrier's solution to include pivots other than the apex. 
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meg, 2.9 Pressure distribution for 5% wedge pivoted 
at apex 


Figure 2.9 presents the component of pressure which is 
out-of-phase with the angle-of-attack for a 5% thick wedge 
oscillating about its apex. The free-stream Mach number is 
mero, Only slightly above that for shock detachment. There- 
fore, the results should highlight possible limitations of 
the potential-flow method. There are two fundamental observa- 


tions. First, the strong effect of thickness, particularly 
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at the low oscillation frequency, and second the performance 
of the present Teipel model, which is considered excellent. 
Calculations on the shape of the oscillating bow shock 
wave are given in fig. 2.10, which shows the dimensionless 
lateral amplitude of the wave at the instant when the aerofoil 
is at maximum incidence. It is clear that the flat-plate 
theory, which forms the basis for Jones' work on the oscil- 
lating biconvex aerofoil [24], fails entirely to predict the 


shock wave locus. 


a theory (T=0) 


Carrier 
(exact 


Present Teipel 





Eon 0 Locus of the oscillating bow wave 
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The main results of the present paragraph show the 
effect of thickness on several stability coefficient 
derivatives for the oscillating wedge. The British notation 


has been used to define the derivatives 


il 
m= In, = f Gee. ear 2. 36a 
0 
1 
m= R, 7 (b-&) eu dé 2.36b 
0 
1 
a (u) 
1 =R, f Ore (Sks 2.36c 
0 


which are, of course, the dynamic and static moment 
derivatives and the lift derivative for the aerofoil, 
respectively. Dynamically unstable oscillations about the 
pivot will thus have -m. < 0 , and the absence OEwampeSieive 
static restoring moment will correspond to “mM. < Oe the 
damping moment derivative —™ is plotted against oscillation 
Frequency in fig. 2.11 for several low supersonic Mach 
numbers. As is well known, the linear theory predicts 
unstable oscillations at low frequency, under certain 
conditions. This is clearly illustrated in the figure for 
each of the Mach numbers considered. Thickness effects are 


seen to increase the possibility of instability for oscillations 
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mide 2,11 Damping moment coefficient for 5% thick 
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of very low frequency. However, as the frequency increases, 
the effect is reversed. At the higher frequencies, in 
accordance with Landahl's observation noted earlier (equ. 2.2), 


there is no significant effect due to thickness. Close to the 
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shock detachment Mach number, thickness effects cause a 
highly nonlinear modulation in the curves of -m, versus iy 
The curves of restoring moment and lift coefficient, shown 
in Rig. 2.12, exhibit similar features. The prediction by 
the present Teipel model of these interesting nonlinearities, 
in close agreement with the exact solution, seems adequately 
to justify the basic isentropic-flow method. Also given in 
figs. 2.11 and 2.12 is Van Dyke's second-order thickness 


solution for very low oscillation frequencies [12]. 
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ieccgse 2. 12 Restoring moment and lift coefficients 
for 5% thick wedge pivoted at apex 


The results discussed above seem to justify Teipel's 


assumption of isentropic flow for the oscillating aerofoil, 


even at Mach numbers close to shock detachment. They also 
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give a useful indication of general thickness effects for the 
wedge oscillating about its apex. Of more practical interest, 
however, are the pressure distributions given in fig. 2.13 
for a 2 1/2% thick biconvex aerofoil pivoted at the leading 


edge, midchord, and trailing edge positions, respectively. 
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As in the case of the wedge, thickness effects are highly 
Significant at the lower of the two frequencies considered. 

At the higher oscillation frequency, however, they seem to be 
important only in bringing about a marked change in the local 
distribution of pressure. The total out-of-phase aerodynamic 
force and moment are evidently in reasonable agreement with 
linear theory. This is shown more clearly in fig. 2.14, which 
illustrates the damping moment coefficient for these three 
torsional modes. The possibility is also indicated in fig. 
2.14 of a significant increase in the stability of the biconvex 
aerofoil, compared with the flat plate, for oscillations about 


rearward pivots. 





fig. 2.14 Damping moment coefficient for 2.5% biconvex aerofoil 
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When the thickness ratios of the wedge and circular 
arc or biconvex aerofoils are 4:1, they will cause the same 
deflection of the airstream at their leading edges and 
consequently they will have the same shock detachment Mach 
number. The damping moment coefficients are therefore 
compared in fig. 2.15 for a 5% thick wedge and a 1 1/4% thick 
biconvex aerofoil. The strong instability of the wedge, 
oscillating slowly about a forward pivot, is shown most 
clearly in this figure. Thickness effects are much less 


important for the ‘equivalent’ biconvex section. 


2.3 Aerofoils in Cascade 

The transonic small perturbation equation is assumed to 
describe the flow throughout the cascade. Consequently, the 
same basic method of solution outlined for the single oscil- 
lating aerofoil will apply. However, in determining the 
fis teady airforces acting on the second blade, there must be 
contended with the fact that the flow immediately upstream 
of the second oscillating bow shock wave is no longer uniform, 
aS is the case for the third blade, and so on. This, the 
remaining problem of significance, is treated in paragraph 
2.3.2. First, however, there must be considered a steady 
oa problem. For the cascade of wedges this is a simple one. 


Yee el, The steady b6 problem 


Consider the region shown in fig. 2.16, downstream 
of the second right-running bow shock wave but upstream of 


the first trailing edge expansion wave. Denote this region 
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fog. 2.16 Preinterference and blade passage zones 


of the blade passage, zone B. The compatibility equation 


for the right-running gB wave through uae is 


Peet) WCE 2. 37a 


a 
Nf We 


while that for the left-running a wave reads 


3/2 


_ 3 8 2 
es pee) ie 0 1 + 5 (VEL) 6 DE STN. 


which follow from the elementary characteristics formulas 
given previously. The state variable xX, which determines 


the steady horizontal velocity Uo, according to equ. 2.12a is 


thus 


3) i 





2s 


ds = ja? = 1) 97? - 3 (ya) 5) 2 oe 


in the preinterference region (zone A), and 


2/3 
7 G ~3)?/4 - cs een) 1 5| 2.38b 


b 
in the blade passage region (zone B). 

Tt 1S worthwhile at this juncture to briefly compare 
this transonic small perturbation theory with the theories 


of Ackeret and Prandtl-Meyer. The pressure coefficient 


Cp = mou. , according to the transonic perturbation theory is 
273 
oS 2 5 | fo?” - = (y+1) 6 - (M | 
a (y+1)M 
2 youd 


on the upper surface of the isolated blade, and 


273 
Cc, = i | [22 = 2 aol) us| ar | 
b Csr 1) MM 
2595 
on the lower surface of the second blade. By taking a 


Binomial theorem expansion for large M (or small 6) there is 


obtained the Ackeret theory, Cc, = DE PMC) and Cc, = 3C, 
a b 
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The transonic and Prandtl-Meyer pressure coefficients, which 
are nonlinear functions of the streamline deflection angle 64, 


are compared with the Ackeret theory below 

















Transonic Small ne - Prandtl 

Pertunmbatton ia Meyer 

C, 0.0486 O204654 0.05065 
a 

Oe Leas 0.13962 O2LGO6G6 


Pea eopmcem woagc Ang leno = Le 


The slopes of the a and 8 characteristics in zone B 
are * , respectively. Also, following the procedure 
for averaging NEw upstream and downstream of the shock, the 
orientation of the right-running bow wave separating this 
mone from zone A is -2/ (fr, + py a 

This essentially completes the solution of the $¢, 
problem. Zone C immediately above the second blade need not 
be considered, since a streamline of the steady flow passing 
through the second left-running bow shock wave clearly 
experiences no further disturbance. Indeed, the upper bow 
wave from the second blade is now an a characteristic of the 
zone A field. Blade-to-blade periodicity of the %6 solution 
in zone C must follow. A typical characteristic mesh with 
the Kantrowitz condition (discussed later) imposed on the 


meeaay anlet flow, is shown in fig. 2.1/7. 


2) 





relations derived previously apply everywhere. 


Since dA/dx is zero. 


ZONE (©) dAc=0-44 


= Wy 
SRN 
Me oN — 


% } Apert! 








fig. 2.17 Characteristic mesh for interacting wedge flows 


2.3.2 The unsteady >1 problem 


The flow 1S assumed to remain supersonic throughout 


tics, respectively 


Z 22 
du, - av, Aj? + i12kM uy/h ax - kM 47% ax = 


ee 2,2 
du, + dv, AJA + i2kM“u,/d dx - k°M$,/A ax 


uy 
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the cascade and consequently the unsteady compatibility 
Indeed, they 
Simplify somewhat for the cascade of oscillating wedges, 


Accordingly, along a ana 8 characteris-— 


= 0 2.40b 





where the state variable X will assume the value ro any) 
zones A and C and Ay in zone B. It is desired to determine 
the change in velocity as the fluid passes through the 
oscillating bow wave from zone A into zone C, assuming 
meemerOpic conditions. Thus, denote u*(x,y,t) and v*(x,y,t) 
the horizontal and vertical velocities immediately upstream 
of P(x,y) on the head wave. The components of u* are L+uy. 
and u,exp(ikt), and for v* they are 6 and v exp (ikt), where 
uy and V, are known functions of x and y from the zone A 
solution of the single aerofoil problem, and 6 is the semi- 
wedge angle. It therefore follows, expressing Uy in terms of 
r, according to equ. 2.12a, that the velocity into the bow 


wave may be written 


Z 
eae tl — 1) ,; 
me =) 1] + + .-- + u, (x,y)e7** 2.4la 
Cy atl.) M 
ee v, (xylen 2.41b 


As with the single biconvex aerofoil (except that the steady 
shock profile is now straight and ee consequently is a 
constant) denote the orientation of the unsteady bow shock 
wave at P(x,y) by Ae as shown in fig. 2.18. Then the 
dimensionless normal and tangential velocity components into 


the shock are 
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mg, 2.18 Geometry of oscillating bow wave for second blade 


< 


: t ae t 
(vtu,)sin(y ty ) (S+v,) cos ly ty ) 2e42a 
— t 7 t 
eS (vtu,) cos (yot+y ) + (S+v,)sin(y.+y ) 2.42b 
Semetting exp(ikt) but keeping in mind that this harmonic 
time dependence applies to all unsteady functions, and where 


the following notation has been used 
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yer eM 


Neglecting second-order products of y' with u anicdueos 


peal! 
equs. 2.42a and 2.42b may be approximated by 


= 
i 


vsiny, = écosy, + (u,siny, - v, cosy, + vy'cosy_) 


2.44a 


= 
lI 


_ fen 
t VCosy, =F ésiny, ! (u, cosy, tv siny, VY siny,) 
2.44b 
The normal component of gas speed Be immediately 
behind the moving bow wave is obtained using the Rankine- 
Hugoniot equation in the form developed previously. Thus, 


employing the present expression for Wr equ. 2.24 becomes 


q as a _ q 
. da ee oe COSY,~V}COSY, dcOSY, Vv 


(y#1) a“U’siny, MSEIa 


Zi 


where v' is the normal velocity of the head wave in stationary 


coordinates. The Binomial theorem yields 
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Z 2 
a 7 2/M _ 2/M 
a = 1 vsiny, a ol 1 >= veosy, ae u 5 : 
s (yh) s y els 
; Lym? 
a | v' 
Veal y2 
2 
2/M 
y-l ‘ieee = _ ee 
+ em (u,siny, V1 COSY, Scosy,) 1 SAT ee: 2.46 


after some Simplification, and where, as previously defined, 
ne is the normal inlet Mach number of the steady flow and 
the quantities y' and v' are of small order compared with 
Yo and U. The unsteady velocities u, and v, are known in 
zone A immediately upstream of the second bow wave. Thus, 
with Ww, = (equ. 2.44b), the components of velocity 
immediately behind the second left-running bow wave can be 
determined. These are the required boundary value quantities 
for the solution of the unsteady characteristics problem in 
zone C. They are obtained from We and Wy, uSing a simple 
orthogonal transformation (2.30). Neglecting such high- 


order terms as y'uys Viuys and y'é, the following simultaneous 


ordinary differential equations can thus be determined 


&> 
_ 
| 


= m, dg/dy + imjg + m,u, + m4Vy 2.47a 


<> 
j 


= n,dg/dy + ingg + n3u, + nyvy 2.47b 
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in terms of the complex lateral amplitude of the oscillating 


bow shock wave g(y). 
eae known from the solution of the 


while the coefficients m and n are 














The velocities u, (x,y) and Vv, (x,y) 


unsteady >) problem, 


known from the solution 


of the steady 6 problem. For the wedge, in which case the 
locus of the steady bow shock wave is a Straight line, the 
mand n are constant. They are given, with ic wo) Bia! 
the general case, as 
m, = BE aii 2 sin? 
IE vtn Yo Yo 
L/M..? 
m, = OS 1 + Sil 
2 y+1 yt cance 
2 
mz = cos“ + va Sin ie 1 - sth: 5 
a (y-1) v 
Sin 2 Sin 2 2/M A 
m = portal) _ yo : 1 - n 
4 ae iL 2 
: y : EY 
5 2.48 
n, = we COS eZ + on sin? 
Ip ve al Yo yt Yo 
Lym? 
n, = as act GOLLY. Sine ¥ 
Z vEL ye O O 
Nn, = mM, 
2 
2 
ie Ssesinay. + ae cos“ l - as 
4 O page fe) i RAS 
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The above equations reduce to those of paragraph 2.1.2 when 

v = 1 and the unsteady perturbation velocities are u, = 0 

and Vi E QO. The unsteady shock locus g(y) and the required 

wy (X,Y) and vV, (x,y) are determined simultaneously in a 

manner Similar to that described previously (by employing 

equ. 2.47 and the a compatibility formula 2.40a). However, 

if the torsional mode of the second blade leads that of the 
first so that its oscillation is described by a> = exp ia (ktth) |; 
where uw is the interblade phase angle, then the unsteady 


flow tangency condition will read 


m7 [ix (o-b) | en’; n= 0 2.49 


which is a simple modification of equ. 2.9b. The vertical 
perturbation velocity v, immediately behind the bow wave 

at the leading edge of the second blade, g(o) = 0, is there- 
fore known and there may be determined consequently the 


corresponding value of uy 


Thus has been presented the basic analysis necessary 
for the determination of the unsteady aerodynamic force on 


the upper surface of the second blade in zone C. Nothing 
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new is involved in extending the procedure to succeeding 
blades. It 1s necessary, however, to indicate how the 
analysis should be modified to yield the appropriate unsteady 
boundary conditions behind the lower bow wave, following 
which the zone B field can be determined. As may easily be 
verified, defining g(y) as an upstream perturbation, the 
changes necessary reduce simply to reversing the algebraic 
Signs of Mo, My, and ny in the shock geometry equations, 
which in this case must be solved in conjunction with the 
8B compatibility relation (2.40b) to yield the required 
U, (X,Y), V5 (x,y). 
2.3.3 Limiting Cases 

There are several interesting limiting situations 
which arise in the case of vanishing blade thickness. With 
v = 1 and M. = 1, the shock geometry equations, written for 
conditions immediately behind the upper and lower waves at 


the leading edge of the second blade, are 


aw 


2V 
= a 2 e t 2 Vea 8) “ 2 1 A 
coal Be Ve he | cos a Sin ‘0 | mL psinzy 











y+1 


> 





< 
| 


2u 
_ eae Oe te 
iL 1 rh Maree hs. ey | sin fea Os al or i> inzy, 


where the upper sign is to be taken for the upper wave and 
where ee is the perturbation in the orientation of the upper 


shock (relative to Vet very close to the leading edge; it 
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memconsidered to be a constant. Similarly 1 denotes a 
perturbation to the lower wave. Both are positive for up- 
stream excursions of the unsteady shock relative to its 
steady-state position. The uy and oH are dimensionless per- 
turbations on the blade immediately behind the bow wave 

and the uy and Vy Similar quantities immediately upstream 

of this point; the latter are zero, of course, for the first 
blade. Now vy is known from the unsteady tangency condition; 
for small amplitude in-phase torsional oscillations about 


the pivot b it is simply 
= -~(l-ikb)a IE ays 


Case (1) Isolated blade 


It follows immediately, with Uae Vy Ot eatha & 





2 
t yl 7 
Y aoa Z ( «sec Vee! ikb) 2293 
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Me th 0 tany, (1-ikb) Zeon 


The significance of these equations is as follows: The 
general shock polar [Ref. 22 page 177], upon assuming the 
turning angle of the flow 6 is so small that there may be 
Wieeecen the relation tan(&-é) = tan & - gsec“E, may in turn 


be written 
sae 6 tan — = sin’¢ = sin*y, Sy. 


where € is the shock angle (measured from the direction of 
the free stream). Writing (-a) for § and with (ry) 

for &, where Ves is second order, equation 2.53 results. The 
isentropic-flow assumptions are thus seen to be consistent 
with the linearized shock polar. The pressure coefficient 
immediately behind the oscillating shock wave at the leading 


edge of the blade is 


C. = + 20/8 (1-1kb) C256 


* 
a result which, it is important to note, is exact in frequency. 


This is simple Ackeret theory with AC,, = 49/8 at maximum 


* 
See equ. 3.13 of the next Section. 
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angle-of-attack, and Be = -4akb/8 when the aerofoil is 
at zero angle-of-attack with the trailing edge moving down. 
Case (a2) Second blade with k = 0 and y = 180° 
The tangency condition for the second blade reduces 
to v = a, while the flow induced by the first blade is 
uy = Oo tany, and ey = -q. Consequently, the shock geometry 


equation (3.5lb) yields 


| 


2 
t — = 
Yi, = a sec oy, (y-1) /2 257A 
oe. Z Z 
‘ioe ac Sec Yo byt3-4cos Vel7s Pee j/h 0. 


which, as shown in fig. 2.20, are positive and negative 
angles, respectively. Substituting equ. 2.57 into the 
shock geometry equation (2.5la) and simplifying results in 


the velocity perturbations -a tany, and 3a tany, on the 


? < 
va \~. 
; N 
P ~ 
; ; ‘, 
\ ne ] Qa. be / 
Ne VN 


70 





upper and lower surfaces respectively of the second blade. 
The pressure loading is therefore — —ieoeay tany.: 
2.3.4 Results 

A mathematical model based on the nonlinear analysis 
has been formulated. During the model validation phase 
several flow solutions were obtained for the oscillating 
wedge cascade for various combinations of the Mach number, 
oscillation frequency, pivot location and thickness-chord 
ratio. These results, however, are of limited practical 
interest until the solution is extended to include the cascade 
exit flow. Consequently, only a representative example of 
these preliminary calculations is given, showing the nonlinear 
effects of thickness on two 5% thick wedges oscillating 
in-phase about a 33% chord pivot at M= 1.2. The inlet 
flow is parallel to the suction surfaces of the wedges; 
otherwise, in the case k = 0, steady equilibrium rotation 
of the compressor blading is not possible. Kantrowitz 
gives a clear explanation of this unique incidence condition 
in Ref. 25. 

The results, for several reduced frequencies, are 
Beevem in fig, 2.21. For the isolated blade, at k = 0.1, 
the strong destabilizing influence of thickness on the 
out-of-phase air loads is clear. For the second blade, 
however, depending on the outcome of calculations downstream 
of € = a, which are beyond the scope of the present inves- 


tigation, it seems highly likely from the results in fig. 2.2l 
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that thickness effects may be stabilizing. The magnitude 
of the nonlinearity in the case of the second blade is 
certainly reduced. At higher oscillation frequencies, 
nonlinear thickness effects are seen to cause only local 


variations in pressure loading. 
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eee over SUPERSONIC CASCADE WITH A SUBSONIC LEADING EDGE. 
- Pee oOo mieORY EN FREQUENCY 


The finite rectilinear cascade was originally solved 
during the course of this investigation to first order in 
oscillation frequency, using both Laplace transformation 
and by generalization of Sauer's solution for the slowly 
oscillating aerofoil in unbounded flow [18]. However, the 
pressure loading was found to remain constant for all the 
blades of the cascade (except the first) while the unsteady 
surface pressures, for the case of oscillations with zero 
phase difference between the blades, diverged in the far 
field. Kurosaka obtains a similar unbounded result for the 
infinite cascade [9]. The work of the present section was 
undertaken mainly to determine if these limitations of the 
first-order analysis could be avoided by seeking a solution 
or higher order in oscillation frequency. 

The velocity potential is first derived exactly to third 
order in frequency for a single aerofoil oscillating in 
unbounded supersonic flow. The cascade with two blades is 
then solved by requiring continuity of the potential along 
the left and right running bow waves emanating from the 
leading edge of the second blade and by suitably modifying 
the flow tangency condition to account for the downwash 
induced by the first blade. The general case of n-blades is 


then given. 
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Be Isolated Blade in Unbounded Flow 


The linearized equations describing the flow field, 


unsteady flow tangency condition and pressure coefficient 


Here che oscillating rectilinear cascade shown in fig. 3.1 


are 
y 
2 376 gee 
9Xx° og? 
ae 
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minere Y defines the aerofoil surface. Consider small 
amplitude torsional oscillations a) = exp(iwT) about the 


pivot x = b and introduce the dimensionless variables 


x 
i 
alx 
Ket 
I 


The above equations, written for the flow in the preinter- 


ference zone of the first blade, become 





2 2 
0° o do do 
p 1 eee 2 1 
B 5 a so) eee sy Be 
OX oy 
9>4 
—— = - [1 + ik(x-b)] ; y= 0 3.3b 
oy 
oo) : 
P = -2 Frage te 1k>)] 3 3C 


The partial differential equation for the unsteady 
velocity potential may be transformed to a more convenient 


form through 


$4 = ye °'™ 3.4a 
2 
ae 
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Maecerms Of the modified potential there results 








2 2 
0 W o° ww 
2 il 2 
eo + 4, - 0 73 
3 oy 
ee ou 7 y = 0 Go. ob 
Cc, = Des 3 - | v5 3.5¢ 
At M 
where the downwash velocity is 
2 
Wy (x) = OR 1+i - (x=-b) 3.6 
M 


Taking the Laplace transform of equ. 3.5a with respect 
to the streamwise coordinate x and noting that the initial 


Bondi tions p, (0,y) and 9p, /ox(0,y) are zero 


av 


— iP) - 87 (p? + 02m?) v1 (p) = 0 3.7 
3Y 


vy (P) denoting the transformed dependent variable. The 


Sorut.on is 
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L 
“By (p +0" /M*) aby (p’+e*/m*)? 


¥, (Pp) = A, (ple 1 B, (p) 


where the condition B, (p) = 0 must be imposed to ensure 
that the flow remains undisturbed upstream of the left- 
running bow wave x-fy = 0. The constant A, (p) follows from 
the flow tangency condition of equ. 3.5b, whereupon the 


modified potential becomes 


w.(p) = ie aby (p?+0*/M*) * 3 3.9 
; B(p* + 8°/M’) * 


Inversion is accomplished using standard tables together 


with the convolution theorem 


S 


f w, (u) J, (ery ru) du Sela 
0 


Py (x,y) 


Rie 


where Jo meomenee Bessel Function of the first kind 


Ja (ry /u) = J = (x-u)i- - B7y- SOD 
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It is now to be assumed that the solution may be expressed 
as a power series based on frequency up to terms involving the 
See Ot the £requency parameter @. Accordingly, expanding the 
mieegrand Of equ. 3.10 and performing the integration 
2 


- 2 
¥) (x,y) = a,s + ids[b stb] + 6°s" [c)stcj+c,x] 


WAS 2 
+ 107s [d,s tdostx(d,std,) ] Seba 

where 

a, = 1/8 

b, = (2M*-1) /2m7g d, = -(8M°-12M+3) /4emg 

b. = -3°b/M“ 2 d, = g*b/6M*B 

c, = -B°/2M8 d, = (2M*-1)/12M"g 

Sy = gb /2M" B dy = p%b/4M"8 

5 Zz 

C2 = -1/4M" 8 Bel? 

and y, (x,y) = 0 for s < 0. The upper surface pressure 


Seerficient for the flat plate oscillating in unbounded 
Supersonic flow, exact to the third power of the oscillation 
jmeeavency, Lollows by substituting equ. 3.l2a into equ. 3.5c 


with k = g76/M2 
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eee Cascade with Two Blades 


Introduce the coordinates & and n, where 


and consider the problem of determining the velocity 
potential o5(&,n) in the preinterference zone immediately 
above the second blade (the shaded region in fig. 3.1). The 


differential equation and tangency condition are 





2 
a“ > a“ > o> 
t—* - —2 = w?m?o, - taku? 54 3.15a 
aE an 
a¢ | 
ar = me NEL + ik(é-b)] 0 n= 0 3.15b 


where the oscillation of the second blade, upon introducing 
the interblade phase angle wu, is described as A. = a,exp(in).- 
The potential must be continuous along the second left- 
running Mach wave separating the preinterference zones of the 


first and second blades. Accordingly 
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$5(¢=0) = 9, (s=m) = py, (s=me*°* 3.16 
Consider the following transformation of dependent 
variable 
=16 
o> (Ern) = p> l(Emnje = ea 
where, as shown in fig.3.l, x = —+l. This is different 
than the transformation of equ. 3.4 on account of the phase 
iaigeexp(—i61), which ensures that along the second bow wave 


@e-s0)(Or S=m) the continuity of ¢ implies the continuity 


of yw. Accordingly, there may be written 
Wo (c=0) = py (s=m) oes 


In terms of the modified potential, the differential equation, 


tangency condition and pressure coefficient are now 


2 2 
0 ow 0 WwW 
2 Z 2 86,2 _ 

ae oo () a) = «CO 3.19a 
3Eo 5 2 M 2 
9V5 | 
7 = —w, (é) exp i(utél) Sg Ae) 
aw 

= 18 _ 2 {. 

Cp, = 38 p> 2 eXp ifo(e+1) +4] 3.19¢ 
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where, by introducing exp(-iu) in equ. 3.19c, the pressure 
coefficient is given at the instant when the second blade 
is at maximum angle-of-attack. 


Take the Laplace transform in the streamwise sense 


ayy (P) De eo 
5 - Bo(p +6" /M")Y., (p) = £(p-n) £210 





dn 


where the non-homogeneous form results from the fact that 
the flow along the line €=0 is now disturbed by the first 
mace, The appropriate initial conditions could, in fact, 
be obtained for d< yn < m/f’ from equ. 3.12a, while for 

n >m/gsp (free-stream) they are zero. It may be noted in 
passing that m, shown in fig. 3.l, 1s the chordwise stagger 
of adjacent leading edge Mach lines in the cascade and that 
the limit m=0 corresponds to the sonic leading-edge locus 
Momartlon. A particular integral for the present problem 
could thus be obtained by employing the Heavyside unit step 
function, as in Kurosaka [9]. In this case an integral may 
be written down immediately from the continuity requirement 
of equ. 2.18 


7 Da: 
Yo (P) = Wz (p) + Ay (pe ON(P FO“) 


aS sa es 
i B, (p) e6N (P Ohm) 


Further, since the inverse of the second term on the right 
vanishes for ¢<0O , there must follow Bo (p) = 0 for 


Wo (G=0) = Y, (s=m) . 
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ime constant A, (p) follows from the downwash condition 
imposed by equ. 3.19b 


~w, (pet (UF0L) 5 aa oe 


Us Dae ies: 
whereupon the general solution of equ. 3.20 is 
L (ut 
we Gee Po a Ufa) ee 
_ 1 SD se 9, 
Wo (Pp) ¥, (Pp) + OEP e 
{p° + 9° /M™) ? 
Seco 
ov, 

where vy" = ae ; & = utl , y=d SZ 


mye inverse is 


C 
Wo (En) = Y) (xy) 4 ~ ff [mer OD ayy J, (Grn, u) du 
0 


3.24 
where the integral vanishes for ¢ < 0 
Recalling equ. 3.10a, this may be written in the form 
L(u+6l1 
Wo (Em) = wy (xy) + vp (Ene™ FEY) 4 r(E,n) 3.25a 


where I(&,n) is to be determined by integrating the product 
of the Bessel function and the downwash function for the 
Eerst blade Vie ViltedtSekhown £Lom the solution of the 


preceding paragraph. That is 
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C 


iCesn) = 2 / W194 (6-1 ,u) du 2A 35) 6 
0 
Neglecting terms of order ae in the expansion of the Bessel 
mone tion 
G 
mee,n) = — Jha + i9[ 2b, (wm) +b, | 


0 


+ ¢? (atm | 36, (utm) + 2c, + 265 (ur) 


+ £0? (utm)| 4a, (atm)? + 3a, (utm) +34 (we) (UF 


+ 2d, (at) | 


du 30 





Performing the integration and defining an interference 


potential Violen), Pears se ounGd sthat 


i(u+81) 


vo(E,n) = ¥, (Ene + W)5(6-N) 5 Bie 


minere the interference term is 
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: Z 2 = 
P12 6Ern) = a,m + iom|by mtb, | = 8 m[cym toymtc, (CTHmx) | 


+ i0°m|dym°+d,m7+dym(3/2 oE4mx) +, (cE +m) | 


Gave + 6h Bz 


Equation 3.27 is the required expression for the modified 
Velocity potential, exact to the third power of the oscil- 
lation frequency parameter 6. It is easily shown to satisfy 
the transformed differential equation and tangency condition, 
while along the leading edge bow wave t = 0 it may be seen 
to reduce to Y1° 

The upper surface pressure coefficient for the second 
blade now follows by operating on by (Ey) according to equ. 
3.19c. The final result, exact to the third power of oscil- 
lation frequency is 


yo 0) 
Cp -- Cp + C e 3.29a 


2 il ie 


where the interference pressure coefficient due to the 
aerodynamic interaction of the first blade may be evaluated 


as 


85 





| 2 
C 2, EGbow hs sue | 6m°E + | mena) +avea + ane? || 


Ie B 28 





SMe | =i AM eye +] -morea ayeMcieaie soe E 


A | =n? or?+4y /3 ~ meatus) - be? (amecas) || 


32295 


being independent of interphase blade angle u and vanishing 
for the cascade with the sonic leading edge locus condition 
(m = 0). The imaginary part of this result, the unsteady 
pressure which is out-of-phase with angle-of-attack, is given 
in fig. 3.2 for two blades oscillating in-phase about the 
pivots at b. The inlet Mach number is 1.1. The results 
labeled characteristics are exact in frequency and were 
derived by neglecting thickness terms in the nonlinear 
oscillating head shock model described in Section 2. 

The low supersonic Mach number strongly emphasizes 
nonlinear frequency effects in these calculations; as shown 
later, the third-order theory gives satisfactory results out 
to much higher frequencies in the higher supersonic range. 


For the cascade considered here, the interference pressure 


coefficient shin compares in magnitude with the pressure 
coefficient eu) for the isolated blade. At k = 0.075, the 
1 
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eg. 3. 2 Upper surface pressure distribution for 
second blade (zero interblade phase angle) 


out-of-phase upper surface lift coefficient for the second 
blade is 1.42; for the same blade oscillating in unbounded 
supersonic flow it is 1.06. 

Turning now to the blade passage zone downstream of the 
second right-running bow wave, for n< 0 and ~ > O (See 
mae oel). The differential equation for the modified 
velocity potential, unsteady tangency condition and pressure 


coefficient given by the set of equs. 3.19 still apply. Also 
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Wo (S=0) = wp, (s=2é+m) 3.30 


ensures continuity of potential along the bow wave 7 = 0. 
In this case the general solution, given by inverting 


equ. 3.21, reduces to 


c 
¥(Em) = V(xy) - of tay ef OM ys Lag ce Few du 
: neu 


where the integral vanishes for € < 0 3; Vy! is given by 

equ. 3.23b; and the over-bar (n) is used to distinguish the 
blade-passage from the preinterference zone. The important 
difference between equs. 3.24 (preinterference zone) and 

3.31 (blade-passage zone) is the Sign attached to the integral; 
however, increased complexity is to be expected in the present 
case because of the reflection of left-running waves from the 
lower surface of the second blade. Some simplification is 
nevertheless possible as only the difference in potential 
through the blade surface is now required. Accordingly 


i (ut+61) 


a (E,0 ) =p, (x,d) - 9, (E,0)e ee) 3.32 


where Vy is given by equ. 3.12a and I(€,0) by equ. 3.26. 


88 





From a previous result 


Pero = tax.) + y,(e,oje 4°) + 1(2,0) 3.258 
whereupon it follows that 
Av, (E) = -2 Fv (oy ef (HOD ; (6,0) | mee 


The pressure loading coefficient is therefore given by 


evaluating 


or = 2p" -- | Ag £{e,0) = 3 | exp -i[ (gen) + 
S234 


the reduction of which, to the expression given below, is 
somewhat simplified by noting that the highest order term 
in I(€,0) occurs as a derivative with respect to €. This 
differentiation is most conveniently performed under the 


integral sign in equ. 3.26. The final result is 
ac, = -2cl4% 4 ac. e 3.35a 


where the interference pressure loading coefficient is 
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AC = get + ao {2 ed - 267m} 
12 1 g 


a Deen 2 2 2 2 2 2 2 ae) 
+ lem E-+{-4m(M"-1) -2M” Bd(2M"-5) ]€+ [-m(2M°-3) +m“Bd+2b8~ ]2M“8d 
28 j 


3.2 
tf (60°42) /36 7+ [2m(4m-3) -848] 67+ [8m(M?—2) -4n Ba (M?-3) 
ZS 
Baia 
+ [-2m*-4n4g*a*/3+4M?mBd (M?—2) ~487a (M7Ba+m) ] Ba} 


Gee SB) 5) 5) 


Considerable reduction is involved in determining this result. 
However, the expression contains only the three basic 
geometric cascade parameters m, d and b (fig. 3.1), together 
with Mach number and frequency. To first order in frequency, 
the preinterference pressure loading of the second of two 
blades oscillating in-phase is thus seen to be independent 
of chordwise coordinate § and pivot location b. This somewhat 
Surprising result applies for low frequency oscillations. To 
the next highest order in frequency, unsteady air forces on 
the second blade depend on both of these quantities. 

A second interesting observation is that when quadratic 
and cubic terms in frequency are neglected, there results, 
for u = 0, simply Sele Bey = 0 and I, (Ap) = 4k/B (M“a/8-m) . 
iors ts) precisely Kurosaka's periodic solution for the pre- 
interference zone of the infinite two-dimensional cascade [9]. 


It is surprising that Kurosaka's solution, the derivation of 
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which requires the inversion of an integral equation, should 
be the same as the solution obtained here for the second 
blade in the finite cascade. An elementary approach to the 
rectilinear periodic cascade problem is given in Section 4. 
Before generalizing the present analysis, fig. 3.3 
summarizes the solution obtained for two blades. The real 
component of generalized aerodynamic pressure acts on the 
blade at maximum angie-of-attack and is therefore in-phase 
with the angle-of-attack. Further, econ) = ~2/8 as the 
frequency of oscillation approaches zero, as given by 
Ackeret theory. The imaginary air forces, however, are 270° 
out-of-phase and therefore act on the blade as it passes 


through zero angle-of-attack with the trailing edge moving down. 






©O CHARACTERISTICS 
~--- FIRST ORDER 
———— THRO. ORDER 


mugee sso Upper and lower surface pressure distributions for 
second blade (zero interblade phase angle, k = 0.2) 


onl 





The range of reduced frequencies over which the present 
solution compares favourably with exact calculations is 
maereated in fig. 3.4. For the cascades shown, the out-of- 
phase pressure loading of the second blade has been integrated 
temeo the point of intersection along the chord of the trailing 
edge expansion wave from the first blade. At M=1.25 the 
third-order theory gives the preinterference lift to a few 
percent for reduced frequencies to 0.4; at M=1.6 the upper 


frequency bound is 0.7. 


@ CHARACTERISTICS 
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—— THIRD ORDER 





micmmeos4) Variation of preinterference lift with reduced 
frequency for second blade (zero interblade phase 


angle) 


a2 





Seo Cascade with n-Blades 


Em 


The analysis is readily extended to the n blade in 


the cascade, With 


meee Preinterference zone of the third blade (fig. 3.1), 
continuity of the potential along the third left-running 
bow wave t = 0, requires 


$4 (1=0) a >> (=m) = > (c=m) e *9* Sow 


The appropriate transformation of » is now 
$3 (ery) a iene 2: 8236 


where x = e+t21. Equations will result which thus have an 
obvious bearing to those of paragraph 3.2. The tangency 
condition (3.19b), for example, becomes in the case of the 
third blade 83/9 = -w, (e)exp2i(y-61) in the plane y = 0. 
Recalling equ. 3.24, the following solution thus applies for 


palery), vValiag for t > 0 


qi 
; OW 
Wa (ery) = > (Een) 2 zf ory (ay oUF OD + Be) |g (ery ew ae 
0 


Seo, 


2) 





But the upwash in the plane of the third blade is now 





from equ. 3.27. The maximum upwash induced in the plane of 
the third blade is thus seen to vary from that induced in 
the plane of the second blade by the amount 9p 5/9n (n=d) , 
which is second order in oscillation frequency. With this 


result, equ. 3.39 may be written 


mel ver(s) 14) : 
W3(ery) = ¥y(En) + ——— ff [oped OD oy, | J (€,¥ 00) 
G 
ic 
a 
+E f syn (n=a) du 3.41 
0 


where, because of the above observation concerning the 
downwash, the expansion of the Bessel function has been 
taken as unity in the last integral. Noting that 

8, >/en (n=d) is independent of streamwise position, and 


employing equs. 3.24 and 3.27 


| MY 
i(u+el) | Wyo (En) + to (n=d) 3.42 


W3lery) = Volerve TT 
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For the preinterference zone above the fourth blade, 


there would be written 


W4(P,r) = palery) + 


\ 
ow . 
- ri [wy (ae! O40) i: x -a)| J, (Pr ru) du 
0 


3.43 


where the meaning of the new coordinates is clear and where 


the upwash in the plane of the fourth blade is 


a (yea) - fe > $22 (nara) = 22 (nea) 3.44 
which may be written 
a a St aa wr (nea) 3.45 
Since 
Ne, 2 Dans eS 
In (n=d) = -mB 7 2056 ania (3nd #244) BaeAG 


from equ. 2.27. It is now easy to show that 


| ay 
Y,(Brr) = ¥y(peret FED) 4 yy (e yn) + AE Ae (nea) 


2t+m . 


where t+T 


See. 





In view of the simple recurrence relation for the upwash 
(equs. 3.40 and 3.45), it is possible to generalize equs. 


3.42 and 3.47; the following result is thereby obtained for 


el 


the preinterference potential of the n blade in the cascade 


355 
on 


HL (are) IL). 


(Eon) = vay (Emde +p,5 (a,b) + [ta-2) c+ | (n=a) 


3.48 


where & and n, with ft = €-8n , have been used as general 


h 


coordinates for the nt blade and Yi92 is to be evaluated at 


a= &+(n-2)1 , b = nt(n-2)d . This simplifies to 


i (ut+é61) 


(6 n) = ¥r-1 (Ene +15 (6-n) 


a2 6° cam* (n-2) { 264na8+(n-1) ml 


“) i6°m?(n-2) { aym/2 [65+ (3n-1) a+ (3n-4) m| 
a a,|26+nag+(n-1) mn] } 3.42%a 


as the modified potential for the nth preimterkrerence zone 


of the finite cascade, with the velocity potential given by 


a Lis) (See (Gel JL) IE) 
(6.n) =, Ceemiie 3.49b 
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It is readily shown that equ. 3.49 satisfies the differential 

equation and boundary conditions described at the beginning 

of this section. The requirement demanding continuity of 

¢(§,n) along the wave ~t=0 may also be shown to be satisfied. 
ele: 


The upper-surface pressure coefficient for the n blade 


is given by evaluating 


(Dye of ule ae . 
Cc = 2}—— WU. —-x~=—] exp -1<4(9[&+(n-1)1) + (n-1)u 
P 2 ere 
ral M 
SO 
where Vn = W660), and cp corresponds to a (max). 
Carrying out the simplification, with k = m“6/8- 
o)  o(w) , ori(n-2)n [of : ze er 
Pp |e P i 
n Ola 1 


where the term of second order in frequency is 


et = k*/8? (n—2) M“m(3m+2d8) 


+ ik?/B’ (n-2) M*m-[m(4M“+1) +3M°a8] & + wh (ENED + m*m? 


- bB2 (3m+2d8) - [m(4M“+1) +2M“dB ] ae} 


3.51b 


which embraces the case n-2. The solution may be written 


24) 





cf _ of) , .(y) -1pH 4 S- alu) .-i(p-l)u 
n 1 dS ES 
p=1 p=3 
Si Royle 
where ey) is the upper-surface pressure coefficient for the 
1 


asolated blade and cfu) 
“1 


coefficient associated with the first and second blades. 


the previously defined interference 


Equation 3.51 is compared with the method of characteristics 


in fig. 3.5. Comments on the solution are given in paragraph 





3.4 below. a 
Third of three blades fo) 
2-0 , aan 
rs S 
Clos | /k= 0:10 
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Fifth of five 
blades 





0-2 Reese - 


0-2 0-4 0-6 0-8 10 ¢ 


fig. 3.5 Upper surface pressure distribution for the aoe 


blade (zer interblade phase angle) 
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The potential difference through the plane of the ne 


blade is 


ay 


i(ut+él ILA 
Ap, (8) = Ayp,-1 (é)e (u ae 2€ ree (n=d) 352 
which gives 
AC. = AC +e tune 3.53a 
Pp Pp Pp 
n n-l 
where the second order term is 
= ke 2. ik? 2 2 2 2 2 2 
AC» = 4 2mamM~ + = madM™~ {[-2 (M°+1) €-2M" lnt+m(4M"-1) +2M" dB-2b8° ] 
B B 
Boob 
Equation 3.53a may be written 
n 
eC ne AG_e 1 (e-L)e 3.53c 
P Po Pp 
n p=3 


mimerne the last term is zero for the case of n=2 and where 


AC, is the loading ccefficient for the second blade. 
2 
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8.4, Further Discussion of the Solution 

PmoOwltlen which 1S exact to first-order in oscillation 
frequency is easily extracted from the more general solution. 
Written for the component of pressure which is out-of-phase 


Peeenetne angle-of-attack, this result for the an blade is 


n-l 
(u) 2 (8: 2km 
= = 1g OS | + ey > | eos joy) 
ge 
5 Ia ©: Se Sin py + 4k Md _ cos 
m De 8B B 8 b 


According to the elementary theory, the out-of-phase surface 
pressure thus exhibits a sustained oscillation with blade 


index, except in the case u=0, when the solution diverges. 


cf) 


E 
n 


from blade-to-~blade. However, the pressure loading or total 


In this case, increases by the constant amount 2km/B- 
damping force remains finite and independent of blade index 


meoeall we. Further, ACy is uniform across the chord of the 
n 


blade and is independent of pivot position. 

The failure of first-order theory in the far field of the 
cascade, for the case when the blades oscillate in-phase, was 
the main reason for seeking the third-order solution. This 
solution is further compared with the method of characteristics 
mmerg. 3.60, Giving out-of-phase preinterference lift for 
several reduced frequencies and interblade phase angles. The 


following observations are to be noted 


100 





© CHARACTERISTICS 
e--- FIRST ORDER 
——— THRO ORDER 





1-50 





1-45 e-—e-—— > —- o-___ 9-9 





© CHARACTERISTICS 
eee FIRST ORDER 
—— THIRD ORDER 





fig.3.6 Variation of preinter- 
ference lift with blade index 
for several interblade phase 
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meqe 2.6 Variation Of preinterference lift with blade index 
for several interblade phase angles 
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(a) for zero interblade phase angle, upper-surface air 
loads given by the WO of characteristics fall below the 
first-order solution, which becomes unbounded as the effect 
of more and more blades is considered. The third-order 
solution is much better. 

(b) the exact results for the surface lift (yp = 0) 
continue to oscillate out to large blade index, although 
the first segments only of these oscillations are shown in 
fig. 3.6 (see, however, fig. 4.2). 

(c) for finite u the approximate theory and the method 
of characteristics predict a strong sinusoidal modulation 
of upper-surface lift with blade index. There is also a 
Sinusoidal modulation in total pressure loading, but of 
considerably reduced magnitude. This important practical 
result is accurately predicted by the third-order theory 

* 


in the near field of the cascade. 


(d) the third-order solution diverges in the far field. 


This is most clearly seen by the occurrence directly of 
blade index in the second-order expressions for iia and 
n 
AC, ‘ 
n 


Notwithstanding the above limitation, the simple closed- 
form solution provides a useful means for estimating generalized 


air forces in the finite cascade over a range of practical 


* e 
The behavior of the exact solution in the increasingly 
far field of the cascade is considered in paragraph 4.2. 


TOZ 





oscillation frequencies. Possible applications include 

(1) preliminary blade design, (11) validation of computer 
models based on exact solutions of the linearized equations 
and, less obviously, (ii1) as a basis for comparison with 

a possible periodic solution for the infinite cascade, 
exact to the third-power of oscillation frequency. In the 
next Section, such a solution is derived which is exact to 
the first-power of oscillation frequency; the surface 
pressure is shown to agree with the mean level of the present 
first-order oscillatory solution. It seems highly likely, 
@eethnis basis, that a periodic solution to the third-power 
in oscillation frequency would agree with the mean level of 


the present divergent oscillatory solution. 


woe 





Meee RLODIC SOLUTION FOR THE UNSTEADY SUPERSONIC CASCADE 
WITH A SUBSONTC LEADING EDGE 


Flutter calculations presently employ generalized air 
forces determined by applying linearized flow theory to 
finite rectilinear cascades [20]. However, no proof has 
been given that these forces are representative of the 
infinite rotor. Indeed, it is frequently very difficult to 
obtain convergence of the surface pressures in practical 
calculations. This has been shown by Verdon [6]. It has 
been repeatedly demonstrated over a wide range of oscilla- 
tion frequencies during the course of the present 
investigation. 

The work of the final Section of this dissertation 
attempts to relate the finite cascade theory of the previous 
Section with an elementary periodic solution derived for 
the infinite cascade. The periodic analysis is based on a 
generalization of Sauers solution for the flat plate oscil- 
lating slowly in unbounded supersonic flow [18]. It is 
exact to the first power of the oscillation frequency. 
Kurosaka gives a difficult mathematical treatment of this 
problem [9]. The present approach, however, is much simpler. 
eee The Periodic Solution 

The dimensionless linearized perturbation equations 
describing the supersonic flow, tangency condition and 


pressure coefficient for the single flat-plate aerofoil are 
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2 2 
0 o> od ao 
2 1 3. 22 2 al 
s) —a = |e a i a2) ee 4.la 
P< a 1 ax 
d¢ 
— Sel Dae ths oie 18))) 1 AE ee 4.1b 
My 
9>, 
“Ph = -2 ae as 1Kk¢) 4.1lc 


aS previously given. The coordinate x is measured downstream 
from the leading edge, b is the pivot and k the reduced 
frequency of oscillation. Upon expressing d, asa first- 
order power series based on frequency, Sauer gives the 


following particular solution 


$, (x,y) = £,(s) + kfhy(s) - iM*y/B£, (s)] 4.2a 
£,(s) = h, (s) = 0: s < Q A2D 
S= x- gy 4.3 


which thus vanishes upstream of the left running leading 
edge bow wave. Apply the flow tangency condition of equ. 


4.1b 


[1+ik(s-b)] = gf! + k (Bhi +im*/B£,) 4.4 
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where the prime denotes differentiation with respect to 
argument and where, Since y = 0, s has been written for x. 


It follows that 
£,(s) = s/B 4.5a 
h,(s) = - is/g>(s/2 + bp“) 4.5b 
The velocity potential for s > 0 is then 
>, (*-Y) = 5 - i my (= + be- + M“gy) 4.6 


and by substituting directly into equ. 4.lc, the upper 
surface pressure coefficient, exact to the first power of 
oscillation frequency 
¢ 4) a ee + i = UDO) RO 4.7 
an B g 
in agreement with the result derived in paragraph 3.1 
using the Laplace transform. 
Consider now the infinite array and denote by x,y the 


* 
coordinates of the flow referenced to the nth blade. At 


* 
No confusion should arise from the use of the same 
x,y coordinates as in the unbounded flow case. 
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time t = 0, as shown in fig. 4.1, the angle-of-attack of 
the at blade is 2s eo that of the ey) 2 blade, 
denoted the gen blade, is O 6 = exp (iy). The preinter- 
ference zone potential o,, (XY) in the shaded region of 


fig. 4.1 clearly must satisfy equ. 4.la. Accordingly 


- re 
De e700) eee ia =) a K| hy (s) iM y/06,,(8)) 4.8 
where s = x - 8y, aS in the previous case. However, while 
C= -Bn=0 m 
. Y 
S=x-By=0 
7 | 


’ 4 ~ W | € 





Picasa, . Cascade geometry 


LO? 





the functions f,(¢s) and h, (s) previously vanished ahead of 


the leading-edge bow wave, this is no longer so. Similarly 


for the preinterference zone immediately above the ae blade 


ey 
; = f eke -iM i 4.9 

Dee n) gh gb n/B (6)| 

6 aCe 1) A ID 
Apply the flow tangency condition to both regions 

my. een 2 

l+ik(s-b) = Bf tk( Bh +1M /Bt.) 4.lla 
lu . we See 

e [a+ix (cb) | = SEE Mis IRI, dil /RE) 4.1l1lb 


in the planes y = 0 (0 < s < 1; S = x) and n = 0 


eC S13; )6CUE 


C), respectively. There immediately 
results for the steady-state and first-order frequency 


terms, respectively 


= i 2? 
f(s) s/B + A, Aa 
z iv, 4.12b 
Ef) t/Be Ag 
h (s) = - 25 (Sipg24m*ga_) - ikB ee 
n Le n n 
5 Sere 2 | Z oe 14 
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where A and B are integration constants. Thus it is 


found that 
. 2 
i S _iks | s Deer 2 ae M 7 
>, (Xr) = Ba, | ie | $4+bg +M Byt+M | TS ya, °, | 
4.13a 
_ Ty Perks Le Tl gy 
4 6Een) E +4 | ce c bR°+M'Bn)e “+M 2a, | 
= bls Me een A 1 3b 
eilaciaa ; 


The requirement for continuity of the velocity potential 


across the bow wave s = mor ¢t = 0 may be written 


which leads to 


by virtue of the steady-state terms in equ. 3.13. 

The velocity potential in the ane preinterference zone 
at the instant when the angle-of-attack of the ae blade is 
a, = 1 has been denoted , (Ss) and at this same time in the 
am zone when Ot = exp(ip) by (0). Thus, at the somewhat 
earlier time when qa. = 1 the potential in the gee zone will 


be @, (or exp(-i1) by virtue of the assumed harmonic time 


Oo 





Seepencgence. For periodicity of the flow 
o at itl 
o,(s) = 6, (6) e 4.16 


Applying this result to equ. 4.13 


-iu 


A, = Se e ay 
and using equ. 4.15 
ee ee 4.18 


ay eal 


The preinterference pressure coefficient for the upper 


th 


surface of the n blade in the infinite cascade, exact to 


the first power of oscillation frequency, follows without 





BH 
cy) = 2 | 2) app | 4.192 
Ph B B a Peal 

which may be written as a modification to C a the upper 


iL 
surface pressure coefficient for the isolated blade in 


unbounded supersonic flow 


: 5 
é = C (u) , i2km/B™ 4.19b 


Pn Pi et -1 
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Consider the flow in the blade passage bounded by the 


right-running bow wave from the oe blade and its reflection 


from the ae blade. The potential may be constructed by 
superimposing upon the previous solution of left-running 


waves a Similar family of right-running waves 
$4 (En) = > (x,y) + Go) “PIs [xg (@)+isPny89g (| 4.20 
c= & + Bn 4.21 
It is seen immediately that 
Sh MON =) 4.22 


continuity of the steady-state potential requiring ve = 9 
when t = 0. With the observation s = ttm, y = nt+d, the 
potential o, (X-y) may be expressed in terms of the local 


coordinates €&,n of the go blade; equ. 4.20 is then 


o 3 Suu = o, (6+) a = - i + [c(m+aan?) +nnen?+c] 


Z 
+ 5, (%) a is cr ~— ng) | 4.23 


miaere C is a constant involving Ae cascade geometry, and 


Mach number. The flow tangency condition applied to the 


h 


lower surface of the qt blade may therefore be written 


sia eet 








- _ 2 
be [rsix(z-»)] ee Beto) +P ik(—S - m) 
M2 
‘Per +k + i 
Bry (BY, sar te 4.24 
where z has been written for —, Since yn = 0. From the 


Steaay-state terms, together with equ. 4.22 


= ean es 
= 4.25 
Tq fs) B 
while from the first-order frequency dependence 
ee: _ A 
i =e togee a: fase (T+bB7) + Mag - nt” | is 
eh 


Hence, the lower surface pressure coefficient for the gq 


blade, from egqus. 4.lc and 4.23 


€ (1) C a 2 fogtkegtik g +m/g- (m+d gM) /8> \ 4.27a 
Pg By q q gq 
(u) in, (a) 4k Ma 2km 
= € + (l-e M)c +i —(——-m +i 
Ph Py BB 3 
4.27b 
, (u) (u) 
where, as previously defined, C and C are the upper 
n the 


surface pressure coefficients for the n blade and for the 
isolated blade, respectively, Finally, there may be written 
for the preinterference pressure loading for the a or 


iggy =e blade, the following expression 


de 





AC =c 1) _ «© () ,ib 4.28a 
Poel Poti Ph 


) 


emg the periodicity of = , and which reduces to 


| 2 
Gl_cae eed i — oa! 4.28b 


These should now be rewritten to correspond to the time 
of maximum angle-of-attack for the (nt1) “8 blade. Thus, 
multiply equ. 4.28b by exp(-iu) and there may be written 
the following general expressions for the preinterference 


pressures in the infinite cascade 


3 
co (4) ae (u) + i 2km/ Bo 4.29 
p aL ea al 
(u) iu 4k M“a fi 
AC 2c (e ite eam) e a a 30 
p Py B B 


4.2 Comparison with Finite Cascade Theor 





The analysis of Section 3 gives the following solution 
for preinterference pressures when terms of quadratic and 


cubic frequency dependence are neglected 


n-l 
ik = oo a8 a > Dp gave PBI 
n it B p=1 
(Ga) 4k M“d -i 
AC 2c ie TA) hg ee ee 
Py BB 


lis 





th as 
mom the n blade of the finite supersonic cascade, with 
u ; 
on ) denoting the upper surface pressure coefficient for the 
isolated blade. Consider equ. 4.31 and for yu # 27m (where 


mis integer) write 


n-1 
( 2k 
let - of] = 28 SE os 


n 


2 _ ; Fae 
_ 2km mL A Sin(n-¥%4s) u 4. 33a 


_ 2km Sin wu Ecos (n=4) 0 
g 3 ZS COSeme 2ucinep/ 2 0 ote 


where it may be noted that interblade phasing introduces into 
the real part of equ. 4.33 terms of first order in frequency. 
The upper surface pressures in the unsteady finite 
cascade are thus predicted on the basis of the elementary 
theory to oscillate indefinitely with blade index, except 
for the case of zero interblade phase angle, when the 
component of unsteady aerodynamic force which is out-of-phase 
with the angle-of-attack evidently diverges. This divergence 
is due to the approximate representation of the velocity 


potential by terms of first order in oscillation frequency, 
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or to some other limitation of the elementary theory, as 
illustrated in fig. 4.2 giving the results of exact calcula- 
tions of unsteady preinterference lift. Considered, with 
zero inter-blade phase angle, are torsional oscillations at 
reduced frequencies of 0.3 and 1.0. The results show an 
almost undamped oscillation in upper and lower surface lift 


Sememnince index out to the 35°" blade. 
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fig. 4.2 Variation of preinterference lift with blade index. 
Method of characteristics with zero interblade 
phase angle 
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On the basis of the low frequency theory of the finite 
cascade, together with the results of fig. 4.2 (see also 
fig. 3.6) it is possible to conclude that surface pressures 
in finite cascade calculations will exhibit a very lightly 
damped oscillation with blade index for most frequencies of 
practical interest in flutter analysis and for all inter- 
blade phase angles. 

Sustained oscillations in surface pressure with blade 
index are not possible in the infinite rotor. Without proof, 
Verdon suggests that a periodic solution may be obtained by 
estimating the mean level of these oscillations [6]. This 
assumption has indeed been the basis for past work on super- 
sonic unstalled compressor blade flutter. As outlined below, 
the present work provides some justification for this approach. 
At low oscillation frequencies, a proof is given that the 
average of finite cascade calculations is in fact the desired 
memrodic solution for the infinite cascade. 

Rewrite equ. 4.29 for the infinite array in terms of real 


and imaginary parts 
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which are precisely the mean values of the two undamped 


oscillatory series given in equ. 4.33 for the finite cascade. 
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The solution for the finite cascade involves the divergent 
series 


nN 


Y } oiPu 


1 


which does not converge to a sum in the ordinary sense of 
the word. However, the series is capable of generalized 
Summation, using the method of arithmetic means [25]. In 


this sense, the Cesa&ro generalized sum is 


> “a Oe > ae 

R, S = 5 + (3 = 4.35a 
HE m=—-© 
y } LH |) Sab oy 

tin . PO ee iu), asia 


where 6 is the delta function. It may thus be said that the 
solution for the finite array converges in the Cesaro sense 
to the desired periodic solution (yu # 0). 

The divergent series for the finite cascade may be 


interpreted conveniently by representing the partial sums 


as points in the complex plane. From equ. 4.33 it is then 


easy to show that these points are equally spaced about the 


iLah 





circumference of a circle which passes through the origin of 
coordinates and whose center is at the point shown (fig. 4.3). 


B@ewamplitude of oscillations is 1/2 cosec 1/2 
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fig. 4.3 Geometrical interpretation of the divergent 
series for the finite cascade array 


The finite and infinite cascade solutions are compared in 


fig. 4.4. Presented against blade index are in-phase and 


out-of-phase components of upper surface lift. The results 
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given include exact calculations based on the method of 
Characteristics which extend to blade index 20. Several 
interblade phase angles are considered. Reduced frequency 
momoe!. Tt is important to note that no periodic solution 
exists according to the present theory when the blades 
oscillate in phase (yw = 0). In fig. 4.4, however, making 
allowance for the expected limitations of first order theory 
meee Lar field of the cascade, it is clear that for u 
greater than some small value, the lift in the finite cascade 
commences to oscillate about the derived periodic solution. 
Further, when the blade oscillations are approximately 
in-phase, the pressure forces in the finite cascade oscillate 
with large period and amplitude. As u approaches 180° the 
period and amplitude become very small. Both the method of 
characteristics and the elementary theory predict this 
behaviour. Note also that the imaginary part of the periodic 
solution is independent of u; however, the real component of 
the pressure varies with u, approaching Ackeret theory as the 
phase angle approaches 180°. 

The unsteady air forces in the finite cascade do not 
therefore converge in the usual sense with increasing blade 
index to the periodic forces predicted by infinite cascade 
theory. Instead, they oscillate indefinitely about the 
periodic solution. Further, while exact calculations indicate 
Eee possibility o£ light damping of these oscillations, the 


rate of convergence is much too slow for practical purposes. 
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In the case of the blade pressure loading, the diffi- 
culties noted above disappear. The AC, GEeeQUss 4.52 SEOrE 
the finite cascade and the AC, eps (2xoGble 3.30 Or Ene an fini te 
cascade awemidentical.,. The fact that AC, is constant for 
all the blades of the cascade (except eee eens was noted 
in Section 3. This finding of the elementary theory is in 


agreement with the results given in fig. 4.5* (see also 


fig. 4.2). 





*Note the large component of lift in these calculations 
which is not associated with the pitching of the aerofoil, 
namely 4 sinu/f 
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5 CONCLUSIONS 


Nonlinear Effects 

A numerical procedure for calculating the nonlinear 
effects of blade thickness in finite supersonic cascades 
with a subsonic leading edge locus is presented. A basic 
feature of the method is the application of the Rankine 
Hugoniot relations to determine unsteady boundary conditions 
immediately downstream of an oscillating bow shock wave. The 
solution is developed downstream of the shock using the two- 
dimensional nonlinear theory of characteristics. For the 
Single oscillating wedge, the solution is in good agreement 
with Carrier's exact solution, even at Mach numbers close to 
shock detachment. There is a significant increase in the 
stability of the biconvex aerofoil, compared with the flat 
plate, for oscillations about a rearward pivot. On the other 
hand, for slow oscillations about a forward pivot, this 
aerofoil exhibits a level of instability which is somewhat 
higher than that given by the classical linear theory. The 
equivalent single wedge, however, is dramatically unstable 
and in this case the pressures show very poor agreement with 
linear theory. The general nonlinear thickness solution for 
the finite cascade reduces correctly in the limit as blade 
thickness vanishes. There is also an indication that thick- 


ness effects may be reduced by cascading. 
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ienear Theory 
The finite flat plate array is solved analytically to 


the third power in oscillation frequency, using the Laplace 
transform. Notwithstanding the divergence of the solution 

in the far field of the cascade, owing to a Bessel function 
expansion which assumes small values of the argument, the 
Simple expressions for pressure distribution have application 
over a useful range of oscillation frequencies. The important 
first-order frequency solution is also derived by generalizing 
Sauer's classical treatment of the oscillating flat plate in 
unbounded supersonic flow. This is compared with a corre- 
sponding result for the infinite eer Te obtained 

using a simple periodic analysis. The pressure loading in 

the two cases is shown to be identical. For the case when 

the blades oscillate with some phase difference, surface 
pressures in the finite cascade exhibit a continuing oscilla- 
tion with blade index. The mean value of this oscillation 


imomonown to be the correct periodic solution. 
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APPENDIX 


Blcernative Derivation of the First-Order Solution 


FOLMENCME Ince Cascade 


A solution for the finite cascade, exact to first-order 
in oscillation frequency, may be obtained by suitably 


extending Sauer's solution for the isolated blade oscillating 





Slowly in unbounded supersonic flow [18]. The flow in the 
aa preinterference zone of the cascade must satisfy 
2 2 
ao a” 3d 
——— Ds k-M~ 6 a i2kM? = A.la 
B yx° 5 Z n OX 


oi (n-l)y 


=e k(x 


where x and y are local coordinates for the nv blade. A 


particular solution is 


o,(x,y) = £,(s) + kth (s) - imMyf,(s)/8] A.2 


where S = x-fy 


and where, for the isolated aerofoil, the functions f(s) 


and h(s) vanish upstream of the left-running bow wave s=0. 


Apply the flow tangency condition 
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Gi ak(sep) jer m7)» 


ee 
Fits 1 ' 
ee cone uh £76) A.4 
where the prime denotes differentiation with respect to 
argument and where, Since y=0, s has been written for x. 


mentollows that 


f(s) = s/per'h7 ln we A.5a 


i (n-1) 


h (s) = eye aye Ne Woy Mpa] +B A.Sb 


where Jin and BL are integration constants. Thus 


3 


¢(xry) = Zen, a - Ik | Geps%an?py) of OD) +WBA, 
g n g 2 n 


6 2 
SiGe 7 6B) A.6 
To determine the pressure coefficient 


, ~i(n-1l 
Cp (ay) = -2[£,'+k(h ‘+i ')) @ fee A.T 


corresponding to maximum angle-of-attack of the ae blade, 
the integration constant AL 1S required. This is obtained 
by matching steady potentials in adjacent preinterference 
regions along the common left-running bow wave. Thus, with 


A,=0, the continuity requirement >, (S=m) = dn+1 (S=0) leads to 
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n-2 
AL, = 5 »» aoe A.8 
p=0 


and the upper surface pressure coefficient for the nth 


blade in the cascade in the preinterference zone is 


n-1 
(u) (u) 2m ~ipu 
on ch + ik S > e A.9 
1 1 es 


where am corresponds to the isolated blade. To O(k) this 
il 

agrees with equ. 3.5lc derived using the Laplace transform. 

If consideration is now given to the flow field bounded 


1216) 


by the right-running bow wave from the (ntl) blade and its 


reflection from the ao blade, then the potential in this 
region may be constructed by superimposing upon the previous 
solution of left-running waves from the a blade an analogous 
family of right waves. The basic procedure 1S given in 
Section 4, where the Sauer approach is applied to the infinite 


cascade. The pressure loading is found to be 


| 2 
Reo = och | oth aa) 4 36k) eee! a.20 
B Py B B 


which is again in agreement with a result of Section 3. 
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